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SOME PRINCIPLES FOR A DATA-PROCESSING 
SYSTEM IN LOGISTICS* 


M. A. Geisler? 
The RAND Corporation 





Some of the principles and concepts of a revised logistics data-processing 
system, proposed for the Air Force and developed at RAND, are dis- 
cussed. These ideas provide a general framework in which a data- 
processing system might be developed that would be able to meet the 
operational needs of the 1960's. The ideas developed in this paper, 
although pointed particularly at the Air Force, might have relevance 
and use for the other military services and for industry in general. 











INTRODUCTION 
Military logistics systems require the use of much paper work and information for 
their performance. Many of the difficulties that have plagued logistics systems in the past 
and up to the present have been the result of weaknesses in the flow of data and information 
throughout the system. These weaknesses have undoubtedly resulted in excessive expenditures 
in the logistics system and consequently a lesser military capability than might otherwise have 
been attained for the amount of resources expended on the logistics system. The rapid changes 
in weapons and therefore in military tactics put an even greater strain on the logistics data- 
processing system so that the assurance of an adequate military capability will require even 
proportionately more expenditures on the logistics system unless adequate data-processing 
systems capable of functioning effectively in the new era can be developed and implemented. 
The RAND Corporation has been doing research on the characteristics and principles 
that might guide the development of a revised data-processing system for the United States 
Air Force in the 1960's, and it might be of general interest to describe some of the results of 
this research in this paper. This data-processing research is part of a comprehensive program 
of logistics research being conducted at RAND, including, as well, supply, maintenance, pro- 
curement, and transportation. The revision of the data-processing system is highly compli- 
cated because it involves the training and education of thousands of people and the introduction 





*Manuscript received June 22, 1956. 

Many of the ideas in this paper have been developed in collaboration with my colleagues at 
RAND. This paper is based on two principal RAND Corporation lications dealing with 
general logistics data processing systems: Research Memo. 1 - I!"Research and Develop- 
ment of a New Data Processing System for Air Foree Logistics," by M. A. Geisler and J. A. 
Postley, and Research Memo. 1417, "A Proposal for a New Air Force Supply Procedure," by 
R. B. McNeill, E. B. Berman, A. J. Clark and H. W. Nelson. 
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of these changes at several thousand locations and organizations. Consequently, the introduc- 
tion of major revisions can occur only during the course of several years, and during this 
transitional period the need to operate two systems creates many more problems than exist 
now or will exist after the new system has been fully introduced. 

In the succeeding sections of this paper I should first like to discuss the operational 
and logistics situation that may exist in the 1960's, for which the new data processing system 
must be designed, the principles of the new data processing system, a brief description of how 
the new system would work, and the plans for testing this system. 


OPERATIONAL AND LOGISTICS SITUATION IN THE 1960's 

The rapid changes in weapon capabilities and design projected for the 1960's must be 
reflected in the logistics system of that period. The growing might of atomic and hydrogen 
weapons is compelling the combat forces to rely heavily upon dispersion and mobility to avoid 
the catastrophe of destruction by a small number of enemy attacks or weapons. 

From a supply standpoint, this means that the demands for parts at each weapon 
location will be different in quantity and character from those at today's large and integrated 
air bases, and it will be difficult to keep adequate stocks of parts at each location. Similarly, 
the present-day depots also will become increasingly vulnerable, and the risks of keeping a 
large fraction of Air Force supply and maintenance resources in a few locations will be very 
great. Thus, the trend may well be to reduce the size of both bases and depots, so as to reduce 
the amount of resources at risk at each location. 

The dispersion of operational units means that the erratic elements found in the present- 
day logistics system will be magnified. The low base-level demands now observed for spare 
parts will be even lower with smaller operating sites [1]. Stocking each of these sites with a 
full range of spare parts would require enormous inventories with little consumption. 

One alternative seems to be a pooling system whereby groups of sites can share in at 
least a portion of stock available at each. This pooling might be arranged in a series of levels 
with certain (high-demand) parts stocked at all sites, other (low-demand) parts stocked one at 
each of a certain number of sites, and so forth. Combined with or alternative to such a system 
might be one in which greater supply reliance is placed on parts in manufacturers’ production 
inventories, with greatly reduced procurement lead times [2]. 

Maintenance policies also will be affected. It will be difficult to sustain continually a 
capability for high-echelon maintenance at each of the dispersed sites. Rather, mobile teams 
may have to move into these sites to make repairs, or the supply system will have to provide 
large numbers of components to the site to replace defective components, which will then be 
shipped elsewhere for repair. The choice between alternatives such as these, which are not 
the only ones, may require major changes in logistics system design. 

There also may be impacts on operational practices. If it is very difficult to maintain 
weapons at dispersed locations, efforts may be made to keep such weapons in a ready condi- 
tion by using them minimally for training and testing purposes, Other weapons located at more 
central locations may be reserved for training purposes, with a rotation of crews between the 
training centers and the dispersed combat-ready sites. 

Weapons systems will need to have a smaller and smaller reaction time to avoid being 
destroyed before they can be used [3]. This would require the logistics system to have 
adequate support in position or to be able to react rapidly in moving initial support to units 
when they deploy and in keeping them supported throughout the campaign. In general, dispersion 
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of stocks and of maintenance capability to a large number of sites, plus the need for fast 
reaction times, will tremendously complicate the inventory control problems, as well as 
distribution and requirements computations and maintenance and transportation scheduling. 

In other words, the conditions of the 1960's will throw a much larger workload on the logistics 
data-processing system. Unless it can handle this more complicated situation, the required 
inventories of parts to make up for the defects in that system will be enormous, and larger 
percentages of aircraft or missiles will tend to be out of commission for parts and for main- 
tenance, and they will be both vulnerable to attack and incapable of offensive action. The data- 
processing system described in this paper is designed to enable the Air Force to cope effec- 
tively with the difficulties of the 1960's. 

The need for the Air Force to respond rapidly to attack means that there will continue 
to be a need for preassembled packages of parts to be prepositioned [4] or moved rapidly to 
the areas of deployment and advanced operations. Even for the relatively stationary organi- 
zations, such as Air Defense units, we cannot expect that the logistics system, with its routine 
resupply, will be functioning in wartime. The size of these preassembled packages and the 
duration of the period they are to support may be different from those used today, but the 
object will be the same as is now sought with the present war reserve system, that of providing 
a relatively inviolate, preassembled, and mobile stock that can support combat organizations 
for an initial period of an emergency, even if the normal peacetime logistics system is dis- 
rupted or lost. 

In summary, then, the logistics system of the ‘sixties will require great dispersion of 
stocks and maintenance capability and great speed of response. Dispersal will, in turn, require 
a data-processing system which can keep accurate and current information on stocks at many 
locations, including the manufacturers' plants. Rapid response will require speed in com- 
pleting logistics computations and instructions to support the fast reaction of the operational 
units. The reduction of logistics-system costs to manageable levels will require better 
logistics policies and more efficient handling of both materiel and information. 


PRINCIPLES OF A NEW DATA-PROCESSING SYSTEM 

An important part of any logistics system revision is to provide a data-processing 
system that can function effectively in the operational environment of the 1960's. All of the 
military services are concerned with this problem, and they are trying to modify their data- 
processing systems appropriately. Undoubtedly, many proposals or solutions to this problem 
will be offered and attempted. There is no obvious or best system for either any one service 
or all services. The subject is a highly complex one involving a mixture of several elements: 
management concepts and rules, communications media, computing equipment, and military 
strategy. There are no straightforward or standard techniques for deciding if one data- 
processing system is better than another system. To a large extent, selection of a particular 
system must be based on general experience, limited development of decision rules and other 
objective guidance, cost considerations, and much qualitative argument. 

This is the process used by RAND in its research toward trying to provide the Air 
Force with an improved data-processing system. The ultimate proof of the system will be in 
its use; all other considerations will leave much uncertainty. Since the Air Force is embark- 
ing, at least in part, on the path described by the principles given below, it was felt that they 
would be of general interest to the readership of a logistics journal. Other publications [5,6] 
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provide much more detail on the characteristics and design of the proposed system than can be 
given in this brief report. The principles describing this system are as follows: 
. Centralized Record Keeping and Data Processing 
. Functional Separation of Stocks and Records 
. Functional Separation of Records and Management 
. Availability of Total Stocks to the System 
. Central Computation of Stock Levels, Requirements, and Other Management 
Calculations 
. Standardization of Input Data from All Organizations to the Data-Processing Centers 
. Differential Processing of Data Applicable to Different Line Items or Organizations 
. Compatibility with Many Different Kinds of Management Systems 
9. Compatibility with Electrical Communication Links 
They will be described in the rest of this section. 


Centralized Record Keeping and Data Processing 

In the revised system there are data-processing centers (DPC's) where policies 
established by the combat command and logistics system can be executed rapidly and unam- 
biguously. These data-processing centers will each handle a large number of line items clas- 
sified into commodity groups such as, airframe, propulsion, armament, guidance, communica- 
tions, general supplies, etc. All the significant information on each item will be available at 
one DPC, These data include, as appropriate, the location, condition, and quantity of the item 
stored at each location in the system until issued to its ultimate user. They also include the 
stock-control level and reorder points for each location and the schedule of due-ins so that the 
DPC can effect routine resupply. Notices of shipment and receipts to the DPC are necessary 
to keep these records current and to provide the data used in many other computations. The 
system-wide asset position of each item can be known currently (rather than after-the -fact) 
only by so centralizing the data on each item. The actual number of DPC's is to be determined 
later. 

The desirability of centralization of information on each item is limited primarily by 
the cost, criticality, the technical limitations of data-processing equipment and communica- 
tions and the need to make the records invulnerable. The electronic state of art is advancing 
so rapidly that the former constraints are continually being relaxed. Various possibilities 
exist for reducing the vulnerability of the centers in general wars, including underground 
location, hardening, duplication of facilities or records, and making wartime logistics 
independent of the centers. 

This last precaution might be accomplished by means of the preassembled packages 
previously discussed or by providing additional inventories at each operating location. The 
latter method may be better suited to stationary organizations, such as Air Defense units. In 
this way, the separate organizations can function independently even if the data-processing 
centers are lost. 


Functional Separation of Stocks and Records 

The fact that the records on each item are centralized implies, of course, that these 
records will be geographically separate from the stocks. This provides flexibility in position- 
ing stocks, since it will be feasible to change the records at the DPC to reflect changes in 
stock locations. Thus stock dispersion and the establishment of mobile forces will be 
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facilitated. Also, it may be feasible to place more stocks at base, which is where demands 
arise, without "losing" these stocks to other Air Force uses. Thus, centralization of records 
can probably be made consistent with decentralization of operations; it may also facilitate it. 
Under the new system, logistics managers may be able to establish detailed policies 
covering many contingencies and have these executed by the DPC's. With the new data- 
processing equipment, it is probable that many transactions that were formerly "exceptions" 
can now become routine, and so the average response time of supply can be reduced. The 
programs prepared for the computer are the translation into machine language of the managers' 
policies and rules. Hence, the DPC's administer the policies of the managers, using the infor- 
mation available on magnetic tape or cards, and thus they can maintain the basic records of the 
logistics system. With this capability, the higher management echelon can be located elsewhere. 


Functional Separation of Records and Management 

Under the present data-processing system, management is dependent to a large extent 
upon the records which can be maintained close at hand. Such records are usually limited in 
scope or detail and many frequently be out of date. Further, since management at each level 
and in each chain of authority requires its own records, there may be a great deal of duplica- 
tion and partial duplication, with possible errors and inconsistencies. 

If in the proposed system there is a single complete set of data on each item con- 
centrated at the DPC, the records will be separated from management as well as from the 
stocks, since many different kinds of managers, including the operational commands, the 
logistics system managers, and other management such as base and depot commanders and 
their staff, will use these data. It is planned that regular reports and specifically requested 
special reports will be made to the different managers in accordance with their needs and 
desires. These reports should be based on one set of up-to-date data, and they should there- 
fore be both timely and mutually consistent. In addition, as will be explained below, the DPC 
should be able to effect policy decisions made by management, even though management is 
removed from the DPC. 

Each organization will, of course, have to keep some records for its own operation; 
for example, stock-location files will have to be maintained at base supply. But there should 
probably be no requirement for one echelon to collect and maintain data which is primarily of 
interest to some other echelon. 


Availability of Total Stocks to the System 
Since it is planned that a DPC will have for each item current knowledge of inventories 
of that stock at all locations, and that therefore it should have the necessary knowledge and 
ability to shift stocks rapidly between locations in response to new needs, these assets can be 
considered as available to the system as a whole, not just to the locations where they happen 
to be stored. This should make possible both more efficient distribution of available stocks 
and reduced requirements [7]. Further, unless stock, wherever located, is available to the 
Air Force as a whole, dispersal may not be economically feasible. Of course, provision must 
also be made for inviolate stocks, such as war reserves, as well as for the incorporation of 
priorities and precedence in the application of the distribution policies. 
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Central Calculation of Stock Levels, Requirements, and Other Management Calculations 

Under the new system, it is planned that managers will establish detailed policies 
covering routine situations, as well as (possibly) many contingencies, so that these can be 
automatically applied by the data-processing centers. This requires the design and specifica- 
tions of management rules and instructions which can be translated into machine language [8]. 
Thus the DPC's should be able to administer the policies of the managers. This should make 
it possible to standardize procedures for many situations that are now considered "excep- 
tional'' and thus make them routine. 

In the distribution of existing stocks in the system, it is planned that normal resupply 
will be initiated by the data-processing Center, based upon its information on each location's 
current stock position compared with its reorder point. Under new distribution rules, the 
reorder point and the reorder quantity also would be computed at the DPC's, with the system 
stockage position and the latest program information taken into account. Bases would thus be 
relieved of preparing routine requisitions. Priority requisitions would be handled in much the 
present manner. 

Central knowledge of world-wide stock positions, consumption experience, and program 
information, as well as faster and improved computing techniques, should permit calculations 
of requirements, maintenance, and procurement schedules that will make for more efficient 
use of available assets and should both increase the effectiveness of logistics support and 
reduce its costs, particularly in procurement and maintenance [9]. 


Standardization of Input Data from the Various Organizations Reporting to the Data- 
Processing Centers 

Uniform procedures and formats for introducing data on issues, receipts, condemnations, 
etc., into the system is necessary to provide greater accuracy and ease of processing. Varied 
reporting systems, with different forms for different transactions, are inevitably accompanied 
by much inaccuracy. The adoption of a ''common language" should also simplify the processing 
of data and thereby reduce its cost. 


Differential Processing of Data Applicable to Different Line Items or Organizations 

Although this revised system requires standardization of input data, it should be able to 
differentiate considerably in the kinds and amounts of processing it gives to different line items 
or to support for different organizations, because the processing work is done at relatively few 
locations and the computers are expected to have much capability and flexibility. This new 
system should also permit changes in management policies to be made comparatively rapidly, 
since such changes need be made only at a few locations by the introduction of a new program 
tape. Consistent treatment of like items or like situations should be facilitated throughout the 
system. Further, because policies must be completely and unambiguously defined in the 
machine code inserted into the computer, greater consistency in the application of the policies 
will be promoted. 


Compatibility with Many Different Kinds of Management 

It should be easier for different kinds of management to work together under the new 
data-processing system. This is made possible by the functional separation of management 
(policy making) and the records and computers at the DPC (administration of policies). Thus, 
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both a weapons-system manager and a property-class manager* can function simultaneously, 
each providing his appropriate guidance to the data-processing centers [10] and, in turn, 
receiving the reports required for his special responsibilities. It is likely that only by this 
functional separation can management concepts dealing with the same kinds of resources, such 
as weapon-system and property-class management, be made compatible. ; : 

Aggregation, for each item, of all significant information at a single point impl. :s that 
virtually any desirable analysis of this information can be made. The DPC is only a service 
organization to the various management entities. Thus, it aggregates data and makes analyses 
at the request of the managements, furnishing each with the reports desired. 

A possible limitation to this is, of course, the DPC's capacity. For example, a 
weapons-system manager may have occasion to make a mathematical analysis of the sensitivity 
of parts demands with respect to, say, the nature of aircraft missions flown. It may be that the 
best way to accomplish this is for the weapons-system manager to make the analysis on his 
own equipment since it is so special-purpose. 


The System will be Compatible with Electrical Communication Links 

High-speed data flow is an important part of a responsive system, but the desirable 
degree of communications speed depends on the cost, error freedom, physical: invulnerability, 
and intelligence security of the means employed. Since high-speed data flow must be accom- 
plished by some electrical methods, it is important to assess the performance, in each of these 
respects, of teletype, radio, and "transceiver" equipment. 

It is proposed that transaction information be introduced at the DPC's in card form 
(where it is transcribed to magnetic tape) and that this information be communicated either by 
card-to-card "transceivers" or courier pouch, depending on urgency and security. It is essen- 
tial that the support system, including the DPC's, be able to operate in limited wars and in 
overseas theaters. Experience indicates that the fewer human beings involved in communica - 
tions, the fewer the errors, and so it is important to mechanize supply communications. It is 
also important that any electrical means of communication include verification devices, so that 
the number of undetected errors is negligible. 


THE NEW SYSTEM 

The principles just described are embodied in a physical and systems arrangement 
sketched very diagrammatically in Figure 1. The core of the revised system is the data 
processing centers, each "prime" for many thousands of items, where the system records are 
kept, the logistics computations made, and the management reports and operation instructions 
prepared. 

The center is envisaged as being a new kind of logistics organization that is functionally 
distinct from present or other new logistics organizations, but linked through appropriate 
communication ties with the operating, management, and command echelons. This revised 
system does not assume that the present logistics organization will be carried forward into 
the next decade without major changes (although it is consistent with much of the present 





*The weapon-system manager is responsible for ensuring that adequate spare parts are available 
to aparticular weapon system, whereas the property-class manager is responsible for ensuring 
that enough of a particular spare part in his property class is available to all weapon systems 
as a whole. The problem of coordination arises when there is not enough in total to meet all 

needs. 
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Figure 1 - Illustrative data-processing system 


structure), and so it is designed to be consistent with the existence of numerous dispersed 
storage and repair sites, perhaps outside the present depot systems. It can also serve a large 
number of small operating bases. It should facilitate greater reliance upon contractor support 
of weapons. It should permit the merging of weapons system and commodity management. 

Three of the most important principles of the new system are indicated in the figure: 
centralization of records and data-processing, separation of stocks and records, and separation 
of management and records. 

The top line indicates several kinds of command and management groups that will direct 
the DPC. Each may provide a different kind of guidance. For example, the combat command 
may state the precedence of its organization. The weapons-system manager may send future 
program data for use by the centers in their calculations. The property-class manager will 
send new catalogue information. Air Force Headquarters may provide budget guidance. Thus 
(principle 8) the system is compatible with many different kinds of Air Force management. 

In return, the management groups will receive a variety of reports on the status of 
inventories for each weapon system or property class, on financial accounting, on weak spots 
in supply, on major changes in consumption patterns, on maintenance and procurement status, 
etc. These reports will all be developed from the same basic data, as is implied by principles 
1 and 6 on the centralization of records and standardization of data. 

The second line of Figure 1 contains the data-processing center. The center should be 
first of all a great store of quickly accessible and currently valid logistics information. It 
should also prepare and disseminate routine operating instructions and management reports. 
Lastly, it should be able to make many of the computations necessary for the administration 
of the logistics system according to policies established by management. 
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The function of the data-processing centers will be to receive, maintain, process, and 
transmit the data needed by the system for inventory control and to compute stock levels and 
reorder points, repair schedules, and the other data described above. 

In addition to the management groups, the DPC has relations with bases, stock points, 
repair facilities, and contractors, as shown in the figure. The relation of the DPC to each of 
these types of organizations will now be discussed. 

The DPC will receive for the line items for which it is responsible notification of 
transactions from each base. With these, it will up-date the balances maintained at the DPC 
for each location and will accumulate historical data to be used in estimating future demand. 
This feature of the system is in accord with principle 1. 

The DPC will also have for each location the reorder points that trigger resupply 
actions, and as base balances fall to these points, automatic resupply will be initiated by the 
DPC for the reorder amount, also contained in its records. Thus, the bases should never have 
to initiate routine resupply actions, although they can initiate requests for priority or special 
needs not met by the routine resupply. The DPC will also compute the reorder points and 
reorder quantities for each location, using the accumulated historical consumption data, future 
program information, and knowledge of the amount of stock available in the system. These 
functions of the DPC carry out principles 4 and 5. 

The DPC will maintain balances not only of the bases but of the stock points (or storage 
sites) shown in Figure 1. It will send operating instructions to these points, such as shipping 
orders calling for the movement of stock to other points in the system or due-in notices to 
inform the stock point of expected shipments for storage there. When a need exists, the DPC 
may also send shipping orders to a base to direct movement of stock to some other point in the 
system. 

The DPC will compute the master repair schedules. It will also receive guidance from 
management as to the repair facilities to be used for overhauling particular items on the 
master repair schedule. Therefore, as it receives notices of reparables generated, it can 
direct their disposition either to some particular facility for repair or to disposal. The DPC 
will also receive notices from the repair facility of repairs accomplished, and it can accord- 
ingly direct movement of this serviceable stock to appropriate locations in the system. 

The DPC will, as necessary, make calculations of procurement requirements that will 
be sent to the proper management groups for review and contractual action. The DPC will 
receive reports from contractors on procurement status and will reflect them in distribution 
calculations, in new requirements calculations, and in directing shipments of newly manu- 
factured products to appropriate locations. 

In all these activities, the DPC functions as a center for using and maintaining records. 
It does not make decisions but only follows the detailed instructions and guidance given to it by 
the appropriate management groups. One of the important functions of the reports prepared by 
the DPC is to inform the management groups on how well the guidance provided to it seems to 
be accomplishing the objectives sought. 


TESTING OF REVISED SYSTEM 
The Air Force has undertaken to test the principles and concepts just described at the 
Oklahoma City Air Materiel Area. At the present time, the procedures and rules needed to 
convert the general concepts into an operating system are being developed to support the B-52 
and the KC-135 aircraft. It is expected that sometime during 1958 the test will begin, with 
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Oklahoma City as the data-processing center. The first plan of this test will be to determine 
what is required to maintain inventory records which are remote from the location of the 
inventory and to resupply these locations automatically by following the rules given to the 
data-processing center. This phase is called inventory control. 

Later phases of the test may involve the testing of new distribution rules that take 
advantage of the centralized and current knowledge of inventory status. Eventually, the test 
may be extended to include calculations of new procurement and scheduling of repair at depots 
or other appropriate repair points. 

Some of the principles described in this paper are also being applied in a currently 
operating data-processing system at Oklahoma City Air Materiel Area for administering a set 
of war reserve supply tables located several hundred miles from Oklahoma City. The storage 
site contains no inventory control records. It relies completely on instructions sent to it by 
transceiver from Oklahoma City. These instructions are prepared mechanically on the IBM 
702 data-processing machine located at Oklahoma City [11]. 


CONC LUSION 

Research indicates that the comprehensive data-processing system described here 
should help much to meet the changing operational requirements of the 1960's. The next few 
years will be required to test this revised system, to develop the detailed procedures required 
for its use, and to extend the system throughout the Air Force if the test proves its merits. 
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PURE STRATEGY SOLUTIONS OF BLOTTO GAMES* 


D. W. Blackett 
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Some necessary conditions for a Blotto game to be solvable by interior 
pure strategies are presented. The mathematics involved is little more 
than an exercise in calculus. Nevertheless, a test of the elementary 
conditions given here exposes as incorrect various solutions to Blotto 
games which appear in the classified literature. 











In a Blotto game [1] two players (A and B) contending on n independent battlefields 
(labeled 1, 2,..., n) must distribute their entire forces (F and G units, respectively) to the 
battlefields before knowing the opposing deployment. The payoff (a numerical measure of the 
gain of A or, equivalently, of the loss of B) on the it battlefield is given by a function 
P; (x,y), depending only on the battlefield and the opposing forces committed to that battlefield 
by A and B. The payoff of the game as a whole is the sum of the payoffs on the individual 
battlefields. 

A pure strategy for A is an allocation of his forces to the battlefields. Such an assign- 
ment may be represented by a vector X = (x), Xo,+++5 X,)» where x, represents the forces A 


n 
commits to the jth battlefield. This vector must have x; 2 0 for all i, and >» x= F. Simi- 


i=1 
larly, a pure strategy for B may be described by a vector Y = Vy erecarens Yn» with all ¥,2 0 


n 
and >»: y= G. The payoff when players A and B use the strategies X and Y is 
i=l 


n 
H (X, Y) - 7 Pi (X;, Y,) . 
i=1 


The particular Blotto games considered here are those which may be solved in pure 
strategies. A game is solved by the optimal pure strategies 
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if H(X, Y) < H(X, Y) < H(&, Y) for all strategies X and Y. In other words, the pure strategies 
X and Y solve the game if, by selecting X, A insures himself a gain of at least H (X, Y), 
while A is limited to this gain if B selects Y. Thus the strategies X and Y are optimal if 

K (X) = H(X, Y) has a maximum at X and L(Y) = H(X, Y) has a minimum at Y. 

Since calculus will be used, the functions P; (x,y) are assumed to be defined and to 
have continuous second partial derivatives for all real values of x and y, with 0 <x < F and 
0<y<G. Those strategies X and Y for which x; > 0 and y;> 0 for i=1,...,n will be 
called interior strategies. 

This paper gives necessary conditions which must be satisfied if interior strategies 
X and Y of a Blotto game are optimal. Henceforth X and Y will denote a pair of optimal 
interior strategies. 


n 
Since L(Y) has a minimum at Y and Zz y,=G for all strategies Y, the function 
i=1 


n-l 
i=] 


must have a minimum when y, = Vy pees Vn = Yn-1 . For these values of y,,...,Y,_;> the 


first partial derivatives of Q vanish. Let L, and Q. denote the partial derivatives of L and 


th 


Q with respect to the k” variable evaluated when Y = Y. Now Q, =L,-L,= 0 for 


k=1,...,n-1. Hence, L,=L,=...=L,. Use Py to represent the first partial derivative 
of Py (x,y) with respect to y evaluated when x = x, and y= Vy: Since 


n 
j=l 


L, = Py. and P, PeeS P. = -B8, where 8 is some constant. The derivative Py represents 
the rate at which additional forces of B at the yth battlefield would increase the gain of A. 


The marginal utility of B's forces at the yeh 


battlefield when X= X and Y= Y is -P,. The 
first necessary condition for X and Y to be optimal is that the marginal utility of B's forces 


is the same for all battlefields when X = X and Y = Y. 
A further necessary condition for L(Y) to be a minimum at Y involves the second 
partial derivatives of Q. Let Qi and Li; be the second partial derivatives of Q and L with 


th 


respect to the i~ and je variables evaluated when Y = Y. For a minimum at Y, the 


(n-1) x (n-1) symmetric matrix 
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Qi . . . Q, n-l Liy+ Lin _ ee Ties 
m=| ° ‘ = Lin Leotlan::: Lin 
Q,-1, 1 Q,-1,n-1 Lon Lon i n-1* Lin 


must be positive semidefinite. If Py is the second partial of P. (x,y) with respect to y 
evaluated when x= X, and y=y,, then L,, = Pj’ and 


‘Al 


tt tr 
» + fe Pp 


tt Lal tt 


P o + « Piiat P, 


From the matrix condition may be derived the following condition: 
' 
For any k of the values P; , the symmetric function formed by summing 


all products of k-1 of the k numbers is nonnegative. 
For k = 2 this becomes 


tt 


tt 
P+ P20 


for any selection of i and j. It follows that P, can be negative for, at most, one of the values 


i=1,...,.m. EP, <0, IP | <|P, | i=1,...,m. 

The inequality P;' > 0 means the first partial of P, (x,y) with respect to y is 
decreasing when x = x and y= y;- In other words, the marginal utility of B's forces is 
decreasing and the gains of B satisfy a law of diminishing returns when x = x and y= y;- 


These conditions, combined with those necessary for K (X) to have a maximum at X, 
give the following theorem. 


THEOREM: 
Hypothesis: X and Y are optimal interior pure strategies of the opposing players of a 
Blotto game with continuously twice-differentiable payoff functions. 
Conclusion: When the players use the strategies X and Y, 
1. The marginal utility of A's (of B's) forces is the same at all battlefields. 
2. The gain of A (of B) satisfies a law of diminishing returns on all except 
possibly one battlefield. The absolute rate of change of marginal utility is smallest on 
the exceptional battlefield. 


REFERENCE 
[1] D. W. Blackett, "Some Blotto Games," Nav. Res. Log. Quart. 1, 55-60 (1954). 
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1, INTRODUCTION 

The allocation or assignment problem is one of the more important problems of man- 
agement science. It may be stated very briefly: We are given a system with a number of 
vacant positions and an equal number of available parts, We know how well each part per- 
forms in each position; we wish to assign the parts to the positions so that system perform- 
ance is optimized. 

Applications range far and wide, from employment to aircraft assignment to naval 
overhaul programs. The computational aspects of the problem have been solved, under the 
assumption that a numerical value is associated with each assignment and that the value of the 
system is given by the sum of the values of the individual assignments [1,2]. The crux of the 
problem, therefore, becomes the finding of values to use for the individual assignments. 

Most of the prior literature has either ignored this problem or failed to reacha 
definitive conclusion. In certain applications it may be possible to determine values in a 
unique natural way; the transportation problem is an example [1]. But in most cases the 
problem is much more complex, and there is no natural utility (like dollar cost in the trans- 
portation case) available. The employment problem is a case in point. Here we have a 
number of vacancies and a number of candidates, and we wish to assign the candidates to the 
vacancies. The problem is of the utmost importance; all organizations but the smallest are 
plagued by it. It has been suggested that, as a first approximation, we classify each candidate 
as either "suitable" or "unsuitable" for a job and solve the problem by using the values 1 and 
0, respectively. This approach is manifestly unrealistic; it is probably worse than letting the 
personnel office struggle along as best it can, using nonnumerical subjective evaluation of the 
candidates and the jobs. 





*Manuscript received September 25, 1956. 
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The technique presented in this paper was originally developed for the problem of 
allocating major electronic equipment to Naval Ships [3,4]. For concreteness, we will refer 
to this application throughout the sequel, except in the last section, where we will return to the 
general case. 

The technique is constituted so as to take advantage of any special internal relationships 
that might exist between the positions or the parts. Thus it will vary somewhat from applica- 
tion to application; and even for a given application, several variations may be possible. In 
essence, though, it is completely general; the same basic ideas apply to all allocation problems 
in which the assignment value is not easily measurable or not measurable at all. 

The naval electronics problem is extremely complex. In order not to obscure the main 
outlines of the technique, we have drastically simplified its statement. A more realistic and 
detailed treatment of that particular problem, of great importance in itself, will be given in a 
subsequent paper [5]. Aspects of procurement and allocation have been covered previously in 
the Quarterly [4]. 

The basic idea contained in this paper is that of Jack W. Smith; the authors are respon- 
sible only for mathematical formulation, elaboration, and technique and, even for these, only 
partly. Much credit is also due John P. Mayberry, Norman Shapiro, and George Suzuki, 

The strictly mathematical part of this paper is contained in Sections 4, 5, and 6; 
Sections 2, 3 and 7 contain more explanatory material than strict mathematical analysis. 


2. OBJECTIVES OF THE TECHNIQUE 

Most problems of operations research are handled in the following way: We are given 
some physical system and a physical utility (dollars, expected kill, etc.) which we would like 
to optimize, subject to the restraints of the system. We set up a mathematical model of the 
system that allows us to express the physical utility (objective function) and the restraints 
mathematically, and then we do our best to find that set of parameters that optimize the 
physical utility. 

This method, although it works very well with many OR problems, fails completely with 
ours. To start with, there is no measurable physical utility that we are trying to optimize. We 
say we are trying to optimize military effectiveness, but this is obviously a subjective idea 
rather than an objective, physical, measurable utility. One can always count his dollars to see 
how much profit he's made; but there's no way of measuring how much "military effectiveness" 
is contributed by some radio transmitter. And even if we had a physical utility, the operating 
conditions of the electronic sets are too complex and involve too many unmeasurable and 
unpredictable factors to be describable at all usefully in mathematical terms. 

The navy handles this problem in the following way: It appoints someone familiar with 
navy objectives and experienced with the equipment and the ships involved, it helps him out 
with a few directives, and it tells him to make up an allocation plan. Although a quantitative 
statement of the objective function and the restraints is out of the question, the responsible 
person can qualitatively weigh all the factors and come up with a "reasonable" or "acceptable" 
plan. In fact, this is how most decisions in the military and business world are reached—the 
responsible person comes up with his decision on the basis of qualitative considerations and 
shrewd guesswork rather than by the use of mathematical optimizations. This process is a 
perfectly accepted and valid one. 

This process works very well when the choice the responsible person must make is a 
clear-cut one with more or less clear-cut implications. Thus if the allocation problems under 
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discussion involved but a few sets and ships, there would be no point in trying to improve 
current methods. Actually, these allocation problems involve hundreds of electronic sets and 
ships; and at this point, the responsible person is no longer able to exercise his judgment 
soundly, because he is overwhelmed by the combinatorial difficulties inherent in problems of 
this size. 

The job of solving large allocation problems really is two jobs: the combinatorial or 
mathematical job and the naval judgment job. Under current methods these two are hopelessly 
mixed together, with the result that probably neither one is being done as well as it could be. 
Our objective is to separate the two: the mathematical or combinatorial problems should be 
solved by mathematicians and computing machines, and decisions involving naval judgment 
should be made by naval personnel. As a result of this process, both jobs will be handled 
more efficiently. 

Note that we do not want to eliminate or replace qualitative naval judgment. In fact, it 
forms the basis of our technique. We are simply out to find a more effective way of using it. 





3. A BRIEF DESCRIPTION OF THE TECHNIQUE 

Roughly speaking, the technique works as follows: The Navy appoints an officer or 
board of officers (henceforth called "the board"') to be responsible for making allocation 
decisions. These decisions will come in response to hypothetical allocation problems of a 
very small size, usually involving no more than two sets of electronic equipment and two 
possible ships to which they could be assigned. Because of the small size of the problems, 
the board will not get entangled in combinatorial difficulties and will be able to reacha 
decision based solely on its judgment, experience, and knowledge, taking into account all fac- 
tors that affect the problem, 

These qualitative decisions of the officers are used in the following way: We assumed 
in Section 1 that each alternative in each of the choices made by the board had a numerical 
value, given by the sum of the values of the individual assignments in it. The choices made by 
the officers therefore yield a set of inequalities involving the coefficients we are seeking. 
These inequalities, together with some additional results on the algebraic make-up of the 
individual coefficients, enable us to estimate the value of the coefficients fairly precisely. 

It is important to note that: 

(a) The value of an assignment is not a function of the assignment only, but also of the 
board that makes the decisions, 

(b) The board is never asked for a numerical estimate, only for qualitative yes-or-no 
decisions, Quantitative values are mathematically deduced from these qualitative decisions, 
assuming only the existence and linearity of the value function. There is nothing surprising in 
this; the process of deducing important numerical information from the mere existence of a 
quantity, without any a priori knowledge of its value, is a commonplace in mathematics. 

We would like to stress the subjectivity of the allocation plan resulting from this 
technique. Our results will be based entirely on the judgment, knowledge, and experience of 
the board; all we have done is to find a technique for obtaining an allocation plan from their 
judgment. If we substitute a different board, we will probably get somewhat different alloca- 
tion plans. There is nothing anomalous in this. All executive organizations, whether business, 
governmental, or military, are based on a chain of command, and it is to be expected that dif- 
ferent people within this chain would make different decisions; this does not negate the validity 
of the command system or the applicability of a decision once it is made. 
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4. THE MATHEMATICAL MODEL — BASIC THEOREMS 

The two classes of objects which we are pairing off are electronic sets, and positions, 
i.e., places on the ships where the sets will be put. We wish to find the value of all possible 
assignments of a set to a position. A 

The electronic sets are classified into models; two sets of the same model are essen- 
tially indistinguishable. We define the void model to consist of no set at all. The positions 
are classified into priority groups, depending on how important they are, and into state groups, 
depending on the model already installed which the new set will replace. The void position is . 
no position at all (conceptually the warehouse). A priority-state group is the intersection of a 


priority group and a state group. 
If a is a set and b a position, we will denote by (a, b) the assignment of ato b. Let X be 











a class of electronic sets and Y a class of positions. An allocation plan T defined on (X, Y) is a 


class of individual assignments (a,b), where a« X and be &, and each a and b occurs at most once 
in T. X and Y will not be fixed in the sequel, but will vary in accordance with the context. (a,b) 
itself may also be considered an allocation plan (defined on ({a}, b})). 

Assumption 1: With each allocation plan T and board B, there is associated a numeri- 

cal value v(T;B). 

Intuitively, v (T;B) can be interpreted as the military effectiveness of T in the opinion 
of B. i ( 

Assumption 2: v (T;B) = a v((a, b); B) . 

(a, b)eT 

From the strictly logical viewpoint, assumption 2 could equally well have been stated 
as a definition. It involves important conceptual assumptions, however, and is used not only in 
theoretical work but in the deduction of numerical values from the board's decisions. It there- 
fore seems more honest to state it as an assumption. 

Denote by f(T) the function that associates with each pair consisting of a model M and 
a priority group P the number of sets of M in positions of P after the plan T is carried out; 
in other words, f(T) specifies the cardinality of the new priority-state groups. Of course, we 
are here counting only those positions on which T is defined. 

Assumption 3: If f (T;) = f(T»), then v (Ty; B) =v (To; B). 


This assumption tells us that the only quantity the Board considers significant can be 
deduced from the number of models of each kind that end up in positions of each kind. Which ( 
set replaced which set is of no interest to the board. Of course T, and To are here assumed 
to be defined on the same set of positions. 








THEOREM 4.1: If a is of the same model as the set installed in b, then 
v ((a,b);B)=0. 


PROOF: Let T be a null plan, i.e., a plan containing no assignments at all. Then, by 
the hypothesis, 


f (T) = f (a, b) 
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Hence, by assumption 3, 


v ((a, b); B) = v(T; B) . (*) 


But since T is null, the sum on the right side of assumption 2 is empty, whence 
v(T;B)=0. (**) 
Combining (*) and (**), we obtain our theorem. 


THEOREM 4.2: If ay and ay are of the same model and by and by of the same 
priority-state group, then 


V (ay, by); B) = v (ag, by); B) . 


PROOF: For i=1 and 2, let T, be defined on ({a,,a},{b,,b)})-and T, = {a,)}. 
Then f (T;) =f (To), and hence it follows from assumption 3 that v (Ty; B)=v (To; B). Since T; 
consists of (a;,b,) only, our theorem follows. 


In consequence of 4.2, we can speak of values v (M, A; B) for triples consisting of a 
model M, a priority-state group A, and aboard B. Formally, we have: 


DEFINITION 4.3: If ae M and beA, then 
v (M, A; B) = v ((a, b); B) . 
(The definition is unique because of 4.2.) 
For a given model M, let S(M) denote the set of those positions that have sets of model 


M already installed (before the implementation of the contemplated allocation plan). S(M) isa 
state group; if P is a priority group, then Pf S(M) is a priority-state group. Define 


(4.4) v (My, P, My; B) = v (My, P 1S (M)); B) . 
THEOREM 4.5: v(M,P, M;B) = 0. 


PROOF: Follows at once from theorems 4.1 and 4.2. 
Denote the void model by 0. 


THEOREM 4.6: v(M,, P,M,;B) = v(M,,P, 0;B) - v (Mg, P,0;B) . 


PROOF: If M, = Mp, then 4.6 follows at once from 4.5. Assume M, # M,. Let bd, 
and by be positions in priority group P. Suppose that by, has a set of model My currently 
installed, and that by has a void currently installed. Let a, and ao be sets of model M, and 
M, , respectively. Let T, be the plan that assigns a, to by and a) tob,. Let Ty be the 
plan that assigns a, to by and ap to by. (Both are defined on ({a,a)}, {b,,b,} ).) Then 
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f(T,) (M,, P)=1=f (To) (My, P) ’ 


f (T,) (My, P) =1-f (T.) (M,, P) ? 


and f (T;) zQs f(T.) for all other pairs. Thus f (T;) =f (To), and it follows from assumption 


3 that 
v(T,;B) = v(T,;B). (*) 
Now 
v(M,, P, My; B) + v (Mp, P, 0; B) = v(M,, PS (Mp); B) + v (My, PAS (0); B) 
= v ((@;, by); B) + v((@o, by); B) 
= v(T ,;B) 
= v(T,;B) 
= v ((a;, by); B) + v ((@o, b;); B) 
= v(M,,PS(0);B) + v (My, PS (M,); B) 
= v(M,, P, 0; B) + v (Mp, P, My; B) 
= v (Mj, P, 0; B) 


Subtracting v(M), P,0;B) from both sides, we obtain (4.6). 


(by 4.4) 
(by 4.3) 
(by assumption 2) 
(by (*)) 
(by assumption 2) 
(by 4.3) 
(by 4.4) 


(by 4.5) 


Intuitively, the expression v (M, P,0;B) gives the value of having a set of model M in 
a position of priority P, rather than the value of a replacement. When P is a priority group, 
we will use the expression v (M, P; B) instead of v (M, P,0;B). These values will be called 
basic values. Theorem 4.6 tells us that, in order to find the values, it is sufficient to find the 


basic values. 
THEOREM 4.7: If b is the void position, then 
v((a,b);B)=0. 
PROOF: If T is a null assignment, then 


f(T) =f (a,b) . 


The result now follows from assumptions 2 and 3. (Actually, (a,b) is a null assignment.) 
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Finally, in order to make use of the decisions of the board, we need: 
Assumption 4: If of two plans, T, and To, a board B prefers T, to T,, then 





v (T,; B) > v(To; B). 
If B is indifferent as to which plan is used, then 
v (Tj; B)~ v(To;B), 


where = denotes approximate equality. 

Those readers who are unhappy about the use of = in a formal assumption may replace 
it by an equality sign. To us, the = seems more reasonable. 

We remark that it would have been possible to state and prove the axioms and theorems 
of this section by using only the terminology of set theory (sets, functions, transformations, 
etc.). Such a procedure would have detracted from the intuitive appeal and added nothing to 
the mathematical rigour. 

Throughout this section, we have used the notation v (T; B) to stress the dependence of 
the numerical value on the board. Since this notation is somewhat clumsy, we will henceforth 
write simply v(T), the dependence on B being understood. B will remain fixed throughout the 
discussion. In the new notation, (4.6) becomes 


(4.8) v (My, P, Mp) = v (Mj, P) - v(Mp,P). 


There is one additional assumption, which will be stated in the following section. 


5. DETERMINATION OF THE RELATIVE VALUES WITHIN A PRIORITY GROUP 

Let us fix attention on a certain priority group P. 

Let b be a position in priority group P that has nothing currently installed. By asking 
the board which model it would most like to install in b, which model would be its second 
choice, and so on, we can obtain a list of all the v(M, P) in decreasing order. (Here P is 
fixed and M ranges over all models.) This is called the goodness order for P. Unfortunately, 
the mere linear order on the functional values v (M, P) does not give us significant numerical 
information about them; in other words, it does not determine their scaling. There is, how- 
ever, a theorem that under certain conditions enables us to obtain a function numerically —to 
within an additive and a multiplicative constant—from a linear order on the differences between 
the functional values. In our case, the additive constant is fixed by the requirement that 
installing a void in b have value 0 (Theorem 4.1). The multiplicative constant will be fixed 
if we specify the value of the highest v (M, P) in the goodness order. This will in general 
depend on P; we will denote it by q(P). Define 





(5.1) r (M, P) = v(M, P)/q(P). 


The r(M,P) are called relative basic values, or just relative values. They are used only as a 
convenience, so that the top member of the goodness order should be 1. 
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THEOREM 5.2:* Let g be a strictly decreasing continuous function on a closed 
interval X. Then g is completely determined by the set of all inequalities of the form 


g (x1) - & (>) {2} g (x) - (x4), 


where x,, Xp, Xg, x,4€X. In other words, if g' is a strict!y decreasing continuous function on X 
distinct from g, then there is a quadruple (x1, Xp, X, X,) for which 


g (xy) - € (Kp) > g (Kg) - (x4) 
g' (x) - g' (X_) < g" (xq) - g' (x4) 


In other words, a linear order on the differences between the functional values determines them 
numerically. The principle of Theorem 5.2 provides the key to solve our problem. It suggests 
that, in order to obtain the basic values quantitatively, we qualitatively compare the differences 
between them; these will be called improvement values. But because we do not yet know the 
value of q(P), we can at most hope to obtain relative values by this method. Furthermore, 
since our problem is discrete (the variable M has a finite domain), whereas theorem 5.2 deals 
with a continuous situation, it is not very surprising that we will be able to obtain these relative 
values only approximately from the qualitative comparisons. 

Let Mj, oeey M.. 1 be all the models, arranged so that 





r(M,P),..., 7M), y,P) 


is in decreasing order (M, +1 8 the void model). This can be done because we already know 
the goodness order. In order to obtain a linear order on the relative improvement values 


(5.3) qi; =e (M,, P)-r (M;, P), 


we must ask the Board to consider all the improvements that can be made, and list in order 
the one it considers most important, the one it considers next most important, and so on. To 
be specific, we ask the board to consider assigning model M; toa position in P that has model 
M; already installed. By (4.4), (4.8), (5.1), and (5.3), this assignment has value qa (P)d, . By 
ranging over all i and j and using assumption 4, we get a linear order on the q (P) qi» which 
yields a linear order on the qd; when q (P) is divided out. 

By means of the goodness order, we can tell at once when a qi is positive and when 
negative. The order on the positive improvement values is called the improvement order. 
Let us denote the largest relative improvement value by d;, the next largest by dy, and so on 
down to’ dy. The improvement order can be written in the form of k-1 inequalities 








*Theorem 5.2 is not new, but we have not found it stated or proved in the literature. A proof 
will be given in a subsequent paper. 
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q >& 
(5.4) ad, > dy 
dy > dy 


It turns out that the inequalities (5.4) are by no means independent and that in fact the entire 
improvement order and goodness order can usually be characterized by a fairly small number 
of inequalities, usually few more than n of them [5]. 

We remarked above that one would need a continuum of models in order for the 
improvement order to determine the relative values precisely. As it is, each of the linear 
inequalities constituting the improvement order can be thought of as restricting the relative 
value vector to lie in some portion of (n+1)-dimensional space bounded by a hyperplane. The 
entire improvement order thus restricts the relative value vector to lie within a certain poly- 
hedron. The smaller this polyhedron is, the better the improvement order determines the 
relative value vector. In practice, the centroid of the polyhedron has been used for the relative 
value vector, but there is really no a priori reason to prefer this point to other points within 
the polyhedron. 

The only information used to determine the relative values is the improvement order. 
Thus, if for two priority groups the improvement order turns out to be the same, then the 
relative value vector for these two priority groups will also be the same. In actual fact, it 
turns out that the improvement order is independent of the priority group. (This also makes 
sense intuitively, because the relative value is in a sense a measure of the goodness or 
suitability of a piece of equipment to a particular mission and should not depend on the 
importance of the position to which it is assigned.) We have been led to: 

Assumption 5: r(M,P) is independent of P. 

As a consequence of Assumption 5, we may write 





r(M, P) =¢(M), 


(g standing for goodness") and conclude from (5.1) and (4.8) that 


(5.5) v (M, P) = q (P) g (M) 
and 
(5.6) v (M,, P, My) = a (P) (g (M,) - g (Ma) . 


q(P) is called the priority rating of P, g(M) the goodness rating of M. We have already 
determined the goodness rating of M. It coincides with the relative value of (M, P) for an 
arbitrary P. It remains only to determine the priority ratings. 








6. DETERMINATION OF PRIORITY RATINGS 

Let P' and P! be two priority groups. Let b! and b’ be positions in P! and PJ, 
respectively, having sets of model mu and mu respectively, currently installed. Let M be an 
arbitrary model, and let a be a set of model M. By (5.6) we have 
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v (a, bY) = q (P) @ (Mm) - g (My) 
(6.1) . 
v a,b’) = q(P’) @ (mM) - g()) . 


We may now confront the board with an allocation problem in which there are only 2 positions, 
namely bi and b! » and only one available set, namely a, and ask it which of the two plans 

fa, b!) and fa, b!) it would prefer. According to its answer, we may deduce from (6.1) and 
Assumption 4 which of 


(6.2) a(P*) (¢ (M) - g(m")) {2} a (Pl) (¢ () - g ()) 


holds, whence we may deduce which of 
(6.3) a )/aP) {2} @ am - ¢ (oe) Agim) - g(a) 
holds. We now draw up a list in size order of all the ratios 

( (M) - ¢ (Mm) /g (M) - gM), 


where M, mi and mM! range over all models for which both sides of (6.2) are positive. 

By asking questions of the above type, and using (6.3), we may find between which two 
items of this list, or near which item, the ratio q (P) /q (pd) lies. Often the list is sufficiently 
dense so that fairly sensitive results can be obtained for any given set of goodness ratings as 
obtained in Section 5. Of course, the "uncertainty" interval surrounding the ratings considered 
in this section depends on the size of the polyhedron considered. 

Once we have obtained all the ratios q (p!)/q (p)), we may obtain the q (p}) and q (p)) 
themselves up to a multiplicative constant (which may of course be ignored). In doing this, it 
is well to fix arbitrarily one of the pl set q (pj) = 1, and then use q (p!)/q (pi) as the value of 
q (p4), in order to avoid the "snowballing" of uncertainties that would result if we combine 
q (Py /q cp) and q (p*) /q (pd) to obtain q (Pp) /q (PJ). Also, in making up the list on which our 
questions are based, it is well to try to stay away from ratios involving small denominators, 
as uncertainties in a small denominator are tremendously blown up in the quotient. 


7. RESERVATIONS AND CONCLUSIONS 

1. Use of the technique presented here should be confined to problems which, like the 
naval electronics problem, involve too many unpredictable or unmeasurable factors or are too 
vaguely defined or complex to be amenable to treatment by the conventional methods of Opera- 
tions Research, i.e., by methods involving physical, objective restraints, and utilities. When 
objective methods are available, they are naturally preferable to subjective methods that 
depend on the opinions of an individual or group, no matter how well informed it may be. It is 
only when we are forced to use subjective methods anyway that the technique described here is 
applicable. 

2. The application of our technique to a particular kind of allocation problem (employ- 
ment, naval electronics, aircraft assignment, etc.) involves the assumption that the model is 
meaningful and valid for that kind of problem. In the final analysis, the validity of the model 
for a given kind of problem must be judged by the military or management authorities 
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themselves. There are, however, mathematical tests which we can use to see if the method 
is at all applicable. The most important of these involves the consistency of the numerical 
results obtained. Before the technique is adopted for use with a particular kind of problem, a 
pilot problem should be run in which several sets of questions are asked (of the same board, 

of course). If widely divergent results are obtained, the model is not applicable to that kind of 
problem. If fairly consistent results are obtained, and if the model is acceptable from other 
viewpoints as well, then it is reasonable to assume that the model is applicable. It is then no 
longer necessary to test for consistency every time a particular allocation is determined, 
because consistency is a logical consequence of the model. Nevertheless, it may be wise to 
throw in one or two more questions than are strictly necessary to obtain the value, just to see 
if the answers check. In the case of naval electronics, limited tests have yielded consistent 
results. It has been pointed out that if linearity (assumption 2) holds, then the answers to the 
questions should be independent of what other assignments have already been made. This also 
can and should be tested. 

It should be noted that no assumptions, implicit or explicit, are made about the "nature 
of psychological judgment,"' except insofar as the routine use of any priority list comes under 
this heading. Psychological or not, all the assumptions are explicitly stated, and the only ways 
in which we can check applicability of the model are by examining for acceptability both the 
assumptions themselves and the results obtained by using them. 

3. There are only three basic assumptions necessary for the application of the kind of 
technique presented here. These are assumptions 1, 2, and 4. Assumptions 1 and 4 tell us that 
associated with each allocation plan and board there is a numerical military value that satisfies 
the intuitive notion of the word "'value.'' Assumption 2 gives us a way of going from small 
allocation problems to big ones, which is of course the crux of our technique. It is important 
to note, however, thai the particular form of Assumption 2—linearity—could equally well be 
replaced by some other formula, and the technique would still be applicable, though possibly 
more complicated computationally. 

Assumption 1 and something like assumption 2 are necessary not only for finding the 
coefficients but also for the application of any method that arrives at allocation plans by linear 
or nonlinear programming. As we remarked in Section 1, standard computational techniques 
for solving the allocation problem require assumptions 1 and 2. As for assumption 4, this 
tells us that we may rely on the decision of the board in arriving at a plan; this assumption is 
certainly justified in the problem we are considering here, because, as stated early in this 
section, we have to rely on subjective decisions anyway. We may therefore conclude that if in 
a given kind of allocation problem linear or nonlinear programming is at all applicable, then 
the kind of technique presented here also should be. 

Assumptions 3 and 5 are not basic; they are merely ways of formalizing the structure 
that the real-life problem in fact has. Most real-life allocation problems probably present 
opportunities for simplifications of this kind. Incidentally, inconsistencies in numerical results 
may sometimes be traceable to ignoring the "fine structure" of a problem in the formulation of 
a "type 3'' assumption. If that occurs, the assumption has to be revised accordingly. On the 
other hand, we are sure that there are some allocation problems in which linear or nonlinear 
programming are not applicable at all. 

4. Even within the framework of assumptions 1 to 5, the computational procedures 
outlined in Sections 5 and 6 are by no means the only ones. We are not even sure that they 
are the most efficient, in the sense that they yield the values as precisely and with as few 
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questions as possible. There is another consideration to be taken into account when asking 
questions. This is the need to phrase the questions as realistically as possible, i.e., to present 
actual situations which the Board can readily picture, so as to make it as easy as possible for 
the Board to consider all relevant factors in arriving at an answer. From this point of view, 
also, it is possible that the procedure in Sections 5 and 6 could be improved. On the other 
hand, the methods of Sections 5 and 6 are easily grasped and lead in a rather straightforward 
manner to numerical ratings. 

In alternative procedures, we might want to determine goodness and priority ratings 
simultaneously, for instance, by exclusive use of questions of the type discussed in Section 6. 
Or else, we might use a completely new type of question; for instance, we could attach realistic 
prices to the models, and use procurement-allocation questions, i.e., questions which ask: 
Given these positions, this money, and these models available for procurement, what would you 
buy and how would you allocate it? The resulting inequalities would yield valuable numerical 
information which could be converted into approximate numerical ratings. 

Some of these various computational systems should be used in checking the appli- 
cability of the model to the problem, as outlined in paragraph 2 of this section. 

5. Another cause for concern is possible variation of the value within the restraints 
set by the decisions we have obtained. Let us consider two sets of assignment values, both 
consistent with the answers we have received. In accordance with Assumptions 1 and 2, they 
will each induce value functions, which we will call v, and Vo: By maximizing vy (T) and 
Vo (T) over all allocation plans T, we obtain allocation plans Ty and Ty that are optimal for 
vy and Vo, respectively. Ty and To may differ considerably without causing undue alarm; 
this is because it is not the plan itself that is important but its military effectiveness. Thus if 


(7.1) 


Vo (T;) x V9 (Tp) ’ 


then we probably need not worry. It is only when there is a pair of value functions not satisfy- 
ing (7.1) that we must try to remedy the situation. In that case we must ask more questions, 
or different kinds of questions, in order to reduce the range of permissible values. In general, 
the more information we have, the less uncertainty there is in the numerical ratings (but we 
must be careful not to ask so many questions that the board ceases to think about the answers). 

If these methods do not help, the usefulness of the technique is seriously impaired; but 
even then it is well to remember that if for two plans T; and Ty, we have 


v(T,) > v(Tg) 


for all value functions v consistent with the answers received, then Ty is probably preferable 
to To. 

On the other hand, it may be possible to obtain plans that are essentially unique (in the 
sense of (7.1)), with numerical information of only the most general and cursory kind. Only 
experience with a particular kind of problem can determine which kind and how many questions 
to ask. 

6. Subjective problems of the kind considered here can have no unique "correct" 
solutions. The technique outlined in this paper is just a way to go about finding acceptable 
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allocation plans, by no means the way. In a particular problem, it may be difficult to evaluate 
this technique; but since the allocation involved usually must be performed in one way or 
another, the technique presented here should, practically speaking, be evaluated in comparison 
with the possible alternatives rather than on its own merits alone. 

Morse [6] has stated that in operations research, the scientist supplies the executive 
with a quantitative basis for the qualitative decision the executive will make. Actually, this 
process is typical of operations research problems on or near the operating level, where it is 
possible to use meaningful physical utilities and restraints to reach conclusions that will be 
used in making decisions on a higher level. The problems with which we are concerned are 
usually several stages removed from the operating level. The solutions must therefore make 
use of the qualitative, subjective decisions made on the lower levels. On the other hand, our 
problems are complicated from the combinatorial-mathematical viewpoint and cannot be solved 
efficiently by the executive without mathematical help. As a result, the process described by 
Morse is here reversed; the executive supplies the scientist and the computing machine with 
qualitative decisions which the scientist and the computing machine use to arrive at quantitative 
conclusions. It remains to be seen whether the same general type of thinking can be applied as 
well to other problems (i.e., not necessarily allocation problems) in which quantitative decisions 
on a high level must be based at least in part on subjective, qualitative decisions on lower 
levels. 
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SOME STUDIES ON SHUTTLE AND ASSEMBLY-LINE PROCESSES* 


M. Pollack 
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Two types of flow processes are described and formulated, a shuttle and 
an assembly-line process. Variations of eachare studied for the general 
N-stage condition. The shuttle process is considered for multiple shut- 
tles per stage and the combination of multiple shuttles per stage with 
storage capacity at the nodes. The assembly-line process is considered 
for multiple machines per stage, both with and without setup times 
included. The mathematical method of recurrence relations is used to 
formulate these relations for the various delay and arrival times of the 
items in each process. It is demonstrated that certain types of shuttle 
and assembly-line processes are completely analogous. All known 
solutions of one process may then be applied to the other. 











INTRODUCTION 
There are many types of processes which may be described generally as the flow of 
discrete objects. Some of these have been called "assembly-line," "queuing," 'waiting-line," 


"bottleneck," or "scheduling processes."' Two such, a "shuttle" and an "assembly-line" 
process, are described and formulated in this paper. The shuttle process arises in cargo- 
handling, the assembly-line process in industrial production. Each is studied for the general 
N-stage process. 

Two phases of problems are encountered in the study of discrete flow processes, those 
of description and those of design. This paper is a contribution to the descriptive phase of 
study. In general, description of these processes is difficult. The recurrence relations 
formulated in this paper attempt to present a simple method of description for the more 
complex types of shuttle and assembly-line processes. These recurrence relations permit 
computational solutions of the arrival and delay times associated with both processes 
investigated, for both deterministic and stochastic types. This study does not treat processes 
where different-size loads are handled in different stages, nor does it treat those processes 
where the items are not processed in order throughout the system. 


I, N-STAGE SHUTTLE PROCESSES 
Previous Studies 

A complete description of the "link-node" model of the cargo-handling system has been 
given elsewhere [1,2,9]. This link-node or shuttle process was evolved from a detailed analy- 
sis of the movement of cargo from a wharf to its final placement in the hold of a ship. 
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The simplest type of shuttle process is one in which there is only one shuttle per stage 
and no storage provision at the nodes except, of course, at the last node. "No storage" refers 
to the requirement that each item be transferred directly from one shuttle to the next. This 
type of two- and three-stage process has been treated by Davis and Weinstock [6]. The 
general N-stage process has defied solution but has been simulated on the SWAC calculator at 
UCLA for stochastic versions [5]. The recurrence relations have been formulated by R. 
Bellman [4] for this simple type of three- and general N-stage shuttle process. The more 
complex types of shuttle processes involving storage at the intermediate nodes and also 
considering multiple shuttles in each stage also have defied general solution. This type of 
complex, three-stage, shuttle process has been simulated for stochastic versions by R. R. 
O'Neill [9]. 


Multiple Shuttles per Stage 

It is desired to determine the recurrence equations for the various delay and arrival 
times of the items and shuttles at each node of a multiple-shuttle-per-stage process with N- 
stages, such as is shown in Figure 1. This process is formulated under the conditions that the 
shuttles do not pass each other. ; 


oF , Vr) Fj , Fiat Fy Pu+t 


oO SHUTTLE 
+H —sSINFINITE STORAGE 


Figure 1 


There are N+ 1 fixed positions or "nodes," P,, Po, one Pua 1 The round-trip path 
between P; and P; “2 is denoted as the jth loop or stage. In the je? stage there are stationed 
mj number of shuttles. The function of the shuttle xX, in the first stage is to convey an item 
(arriving from source O) from P, to Py and deliver it to X,, a shuttle in the second stage. 
xX cannot return to P, until X, has accepted the item at Poy. Upon receiving the item from 
X,, X conveys it to Pg, where it is delivered to Xg. Similarly, X cannot return to Po 
until X has accepted the item at Ps, etc. After Xy has picked up an item from Xy 1 at 
Py Xy transports the item to Pue Dp the terminal, and immediately returns empty to Py: 
The shuttles in the j stage are denoted by the symbol x E; The first subscript, j, denotes 


the stage; the second subscript, Ej, distinguishes a particular shuttle. E; = en m)- The 
subscript E, is used to determine which shuttle (in the j*” stage) transports the i*” item 
according to the following relation: 


a) Ey =1- (m,) («), 
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where i=the i” item, i=1,2,....n, 
m; = number of shuttles in the jt stage, 


k = an integer such that 1 <&, <m, . 


Given a set of items, T, with associated transport times, the following quantities 
describe the set: 


X, jg (i) = the time required for X, to transport the i*® item from P, to P,. ;, 
LE; ,E, j Phot 


7 r] Ej 1 


' 
x, Eg,” = the time required for x, E, to return to P; from Py y» after having 


th 
delivered the i~ item to ‘i 
Ki 1, Ej ot 


The actual pick-up and delivery times are assumed to be zero. Or they may be thought to be 
lumped in the transport and return times. The delays are defined as follows: 


pa th 
qj 2,” = the delay to X E, awaiting the i~ item at P; ‘ 


4 (i) = the delay to X waiting to deliver the i‘ item at P 1° 
iE, E; i 


Finally, the following arrival times are defined: 


Tr (i) = the time at which x; E arrives at P, ready to receive the jth item, 
i, j > j j 

8, g. (1) = the time at which X, ,, arrives at P,, , carrying the it? item. 

i, j “ j+ 


Since m; is known for every j, subscript E; is redundant above. It will be retained only when 
required for clarity, otherwise replaced by a star (*). 


In the j® stage there are my shuttles, XH Xo eee x m’ where the general term is 
Xie Since all the shuttles are initially at Pi at time equal to zero, 


(5) T5«(4) = 0 i< m; . 


The time of arrival of the first shuttle, XH at Pj 4 ‘transporting item 1, is 
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The delay, dj, (1), experienced by Xy waiting for shuttle %-1,1 to arrive at Py is defined 
in general as 


(6) jo) = max [85_1 +(1) - rye0), 0]. 


The time of arrival of the second shuttle, Xj at P ei transporting item 2 is 


8;9(2) = max [5j1), dijo (2) + X2(2)| , 


where max [a,b] denotes the algebraically larger of the two values a and b. The above 
relation for s 9 (2) is’seen to be clear if it is remembered that the shuttles may not pass each 
other. The arrival time of Xo is constrained to be the larger of their independent arrival 
times. In general, 


1) sje) = max [554(1 - 1), dy) + x40] , 


i< my; » 
where, for i= 1, 8;+(0) =0. 
The time of arrival of shuttle X1 back at Py» ready to receive item m, +1, is 
Ti (m, +1)= 85 (1) + X) (1) ° 


This is so because there is always at least one shuttle in the j + ith stage which initially is 
always at P = and therefore the delay to %1 at Pi 1 is zero. If not constrained by Sp 
the independent time of arrival of Xe back to Pj is 


8;9(2) + 4; (2) + Xo (2) . 


The delay 459 (2) to Xo at Py + 1 may be zero or some positive number, depending on whether 
or not there is a shuttle waiting to receive item 2. Therefore, the actual constrained time of 
arrival is in general 


“ rj«() = max [Fy+4 - 1), 8;,4(i - my) + 4 jf - m)) + jel - m,)| ‘ 


i>m,. 


The time of arrival of Xj a second time at Pi+ j» Conveying item m; + 1, is 


8;+(m; + 1) = max [sy+(my), Ty+(m, + 1)+ d;4(m; + 1)+ Xj«(m, + 1)] " 
And, in general, 


sje) = max |sy4 - 1), rya(l) + dya(t) + xyoGt)],, 
Ce een 


(9) 
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This equation is also valid for i < mj, since by Eq. (5), 
T5«(i) = 0 for i <m,- 

The delay to a shuttle X« waiting to deliver an item at P; * is defined as 

(10) 4; +(i) = max [Ti+ 1,+@) - Sj«(i), " , 


where 454) =0,i< M419 since there are m, +1 Shuttles initially positioned at P; +1° 


The general jh stage recurrence relations are then (5), (6), (8), (9) and (10). These 
relations may be used for any stage except the 1st and nth, 


For the 1st stage, d,+(t) = 0 iz1,2,...m. 


For the N™ stage, Ana(i) = 0 i= 1,2,....n. 


Multiple Shuttles per Stage with Storage 
This complex shuttle process, shown in Figure 2, operates in identical fashion as the 
one previously described with the added feature of a storage allowance at each node. 


7" : 7 Fi : Fel AN Fel 


's] SHUTTLE 
+H —SsCSINFINITE STORAGE 
alll FINITE STORAGE 


Figure 2 


The quantities which describe the set of items, T, include (2), (3), (4), and also 
(11) *, = the number of items that can be stored at P; +1° 


The recurrence relations which describe the jt stage of this shuttle process consist 
of (5), (6), (8), (9), and the following relation for the delay, A je) , 


A.,(i) = max[r.., 4(i-0,) - s,,(i), 0], 
(12) ia al 


i > M+ * 
where 


(13) A jf) = 0, i<mj,)+ o;- 
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In (13), 4 j«(1) is equal to zero, since the first M41 items are delivered directly to the await- 
ing shuttles in the j+ yth stage; the next °; number of items can then be deposited in storage 
at Pi, 1 In (12), the delay is zero if a shuttle in the j+ ith stage has removed item i -o 


j 
from node Pj +1 before the shuttle in the = stage has arrived with item i. 


For a process with infinite storage at all nodes, the A y+) delay is always zero. 
(14) 4.) = 0, Se éh.@...8,; Jeee...8. 


Il, N-STAGE ASSEMBLY-LINE PROCESSES 
Previous Studies 

The "assembly-line" processes discussed in this section are an idealization of the 
actual behavior of an industrial assembly line. There are a series of stations or "machines" 
where each item, in order, is subjected to some process or modification. After the process on 
an item is completed, the item proceeds to the next machine. Here, either it is immediately 
accepted into the machine or it must wait if there is a waiting line. This type of two- and 
three-stage assembly-line process has been treated by Johnson [8]. Recurrence relations 
have been formulated by R. Bellman [4] for two-, three-, and general N-stage processes. 

As in the case of the shuttle process, the more general N-stage assembly-line process 
has defied solution. Also, the more complex assembly line which includes multiple machines 
per stage and also restrictions on the length of the waiting lines has defied solution. 


Multiple Machines per Stage 

It is desired to determine the recurrence relations for the various delay and arrival 
times of the items at each node of a multiple-machine-per-stage process with N-stages. This 
process is shown in Figure 3. 


G MACHINE 


Figure 3 


There are N+ 1 fixed positions or nodes Pi, Po, «+: Puy te The interval between P 
and P; +1 is denoted as the j® stage. In the je stage are situated M, number of identical 
machines. The machines in the jh stage are denoted by the symbol W. F.’ Similar to the 

> 


j 
shuttle symbol, the second subscript F; distinguishes a particular machine. 


j 


(15) F)=i- (M;) (k) 
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where i = the it” item, tot &...ce, 
M, = number of machines in the jth stage, 


ete 


k = an integer such that 1 < F; < M;. 


Initially, there is a waiting line of items behind "'gate"' P,. Py at time zero, with no 
delay, lets item 1 into machine Wip item 2 into Wie etc., until all M, machines have items 
inthem. "Gate" Pi lets items out of the 1st stage only in their original order. Even if the 
processing on item 2 is finished before that on 1, item 2 will remain in its machine until 1 has 
passed Pi Immediately item 1 has passed PY item 2 will pass through and also let item 3 
through if its processing is completed, and so on. As soon as item 1 has left Wip Py allows 
item M, + 1 to enter W,,. Only after this item has passed P, does P, permit item M, + 2 
to enter Wie (if it is through processing item 2). After the items pass by Pi they imme- 
diately proceed to the waiting line behind P, for the second stage, and so on. 

Given a set of items T', with associated process times, the following quantity describes 
the set: 


(16) Wy F (i) = the time required to process item i on machine F; in the j* stage. 
5 A 
J 


The waiting line delay is 


(17) Dj.) = the delay to the it” item in the waiting line for the j+ 1" stage 


(behind P ), after being processed on W - 
j+1 , F; 


The time of arrival is defined as 


(18) S, y, (i) = the time of arrival of the ith item at the waiting line for the j+ 1°? 
> 


stage (behind P; is »» after being processed on W, ,,. 
ve 


There are M; machines in the jh stage, designated Wy W32 ous Ww; My where the 
> 


general term is W3+ At P,, the delay to items 1 through M; is zero. 


j? 


(19) Dy - 1, #() = 0 i<M,, 


since all M; machines are initially empty. When item 1 is finished processing in Wy1 it 
passes Pi and also arrives behind P; «4 at time 851 (1). Item M, + 1 arrives behind Py 
ready to enter Wy1 at time Ss; & «(My +1). If S51 (1) > S)-1, «(My + 1), then item M; +1 


will experience a delay 


Dj _1, +(My + 1) = max [Sy - 8) «(Mj + 1), 0}: 
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And in general, 
(20) Dj _3, +4) = max [S+4 - Mj) - 8;_1 *(4), | i> M,. 


The time at which item 1 passes P; and also arrives at the line behind P; +138 
‘ 


8511) = S)_4 4) + Dp_y 1G) + wy). 


' at 


If unconstrained by item 1, item 2 could pass P j 


S_1, 2(2) + Dy_y 92) + wj9(2). 
Its actual time of arrival is then 
Sjq(2) = max (8541), 81,202) + Diy, 2(2) + wyp(2)] . 
And, in general, 
(21) Syo(i) = max [Spo(t - 1), 851, x(t) + Dj_y, of) + Wyld] , 
be && ..08. 
where Sp (i) = 0, i= 1,2,....m. 


Multiple Machines per Stage with Setup Times Included 
The set of items +. defined in (16), is modified to include the setup times. 


(22) Wj F,”) = the time required by machine W; r, (machine F; on ha stage) 


to setup for its next item after having processed item i. 


The following recurrence equations describe this process: 


Dj 1, +) = 0, i< M, 
max [Sy4(t - My) + wix(l - My) - Sy _ 1, 4(4), 0] 
i> M; 
Sj4() = max [Sj+4 - 1), 8)_4 of) + Dj _y oi) + wye0)] ; 
Soh, 3 2-08, 
where Sp) (i) = 0, for i = 1, 2,....m. 


Relation (23) differs from (20) since the setup on machine Wis cannot be accomplished 
until item i has passed P)- Only after the setup interval is W3x ready to receive its next item. 












SOME STUDIES ON SHUTTLE AND ASSEMBLY-LINE PROCESSES 


m. ANALOGY BETWEEN SHUTTLE AND ASSEMBLY-LINE PROCESS 

It can be demonstrated that the equations describing the shuttle process in Section I 
are identical to the equations describing the assembly-line process in Section II. The typical 
* stages of each process are shown below in Figures 4 and 5. 








Sh ia 

_ + 
«a x 
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Figure 4 Figure 5 


For the shuttle process, the equations were formulated to describe the activities of the 
shuttles. For the assembly line, the equations were formulated to describe the behavior of the 
items. It is necessary to formulate for the assembly line the equations describing the activities 
of the machines, and for the shuttle process the equations describing the behavior of the items. 
Comparison can then be made of the appropriate sets of equations. 

In addition to the variables defined in Section I, the following quantities are defined for 
the assembly line: 


(25) R; F (i) = the time at which machine W; Fr, is ready to receive item i for 
es yes 
processing. 


dj F” = the delay to Ww; F; awaiting item i, after it has become ready to 
’ ’ 


receive it. 

A, E (i) = the time at which item i is finished processing on W. , F on the 
™ stage. 

55 r” = the delay to Wj, F; awaiting item i to leave, after it is finished 
processing. 


For the shuttle process the following additional quantities are defined: 


(26) aj Ez,” = the "unconstrained" arrival time of x, E, at P; +1 Carryingitem 1. 


6 = Al we 
i, Ez,” the delay to x E, carrying item i to P; 41 due to being "constrained. 


D; Ez,” = the delay to item i, after delivery into the storage at Pi ey by 
, 


ial X, g, waiting to be picked up by its next shuttle in the j+ 1° stage. 
) ? 
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When these recurrence equations are formulated, both the pair of equations describing 
the shuttles and machines and the pair describing the items in each process are found to be 
completely analogous. The analogous quantities for each process are all shown in Table 1. 


TABLE 1 
Analogous Quantities 





Assembly Line Remarks 





Position or node 





Item 








Stage 





Number of stages 





Shuttle or machine 





A particular shuttle or machine 





Number of shuttles or machines in 
a stage 





7 


Set of items 





4,8, 


w (i) 
iF; 


Transport or process time 





tO 
“5, 


W; p (i) 
iF, 


Return or setup time 





r (i) 
iE; 


(i) 
iF; 


Arrival time of shuttle or machine ready 
to receive item 





d. » (i) 
iE, 


d; . (i) 
i, F; 


Delay to shuttle or machine awaiting 
item 





s (i) 
iE, 


S; p (i) 
iF; 


Arrival time of shuttle or machine ready 
to deliver item 





a (i) 
iE, 


A; F, (i) 


Unconstrained arrival time of shuttle or 
machine ready to deliver item 








6. (i) 
J, Ej 


5,7, 


Delay to shuttle or machine due to being 
constrained when delivering item 





D8,” 





Dj, F, (i) 





Delay to item in storage or waiting line 
awaiting its next shuttle or machine 








The assembly-line process was formulated with no restriction on the size of the waiting 
lines before each machine stage. If a restriction is placed on the size of these waiting lines, 
then this modification corresponds to the finite storage provision in the shuttle process. 












yaiting 


nes, 








SOME STUDIES ON SHUTTLE AND ASSEMBLY-LINE PROCESSES 


Iv. COMBINED SHUTTLE AND ASSEMBLY-LINE PROCESS 
Some flow processes combine features of both shuttle and assembly-line processes. 
Such a three-stage process is shown in Figure 6. The first stage, Py ~ Po, has three 
machines and storage capacity for 10 items at 5,. In the second stage, Po - Ps, a single 
shuttle transports items between the two machine stages. There is no storage provision at Pg. 
The shuttle must therefore wait at P, until a machine on the third stage is empty. The third 
stage has four machines, with infinite storage at tiie terminal, P,. As shown in Section 3, a 
stage of a shuttle process can be easily formulated which behaves in an identical manner as an 
assembly line stage. Stage Py - Po in Figure 6 can then be replaced by a three-shuttle stage 
with storage for 10 items. Also, stage Ps - P, can be replaced by a four-shuttle stage, with 
infinite storage at the terminal. This analogous complex shuttle process is shown in Figure 7. 





Figure 7 


The second-stage, single-shuttle process could have been replaced by a single-machine 
stage, thus forming the completely analogous assembly line process. 

Stage-by-stage transformation of any shuttle process by an assembly-line process (or 
vice versa) can thus always be accomplished. Therefore, computational or analytic solutions 
of one flow process may always be applied to any analogous process. 
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A TRANSPORTATION AND PRODUCTION MODEL* 
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A transportation and production model is given in which costs of trans- 
portation and production are considered. There are no restrictions on 
the number of transport paths, and losses in transportation as well as 
those at terminals are incorporated in the model. Also, products of one 
type may be converted into products of another type at knownterminals. 
There are many applications of interest in the use of this model, where 
the demands for products are known and the desired solution allocates 
production in such a manner that the total cost of maintaining the system 
is minimized. 











INTRODUCTION 

A general problem in the economics of production is how to satisfy a demand with least 
cost. If the cost of moving the product from source to consumer is an appreciable part of the 
total cost and if there are a number of different possible routings, the problem also becomes a 
transportation problem. The transportation and production problem formulated here is for the 
least-cost allocation of production among available factories to meet known demands. The 
formulation assumes known demands and known costs of production and transportation, permits 
transportation and production losses, and considers the production of interdependent products. 

Various methods are available for solving a mathematically formulated minimization 
problem. These include direct gradient methods, analytic schemes in which formulas from the 
calculus are used, and, in the case of problems expressed by linear functions only, linear- 
programming methods. When the cost functions for production and transportation and the loss 
formulas are nonlinear and when the problem includes many variables, the known methods of 
solution become cumbersome to use, and actual solutions often are not feasible. Usually, 
problems of this type are solved by making linear approximations, by lumping terms, by 
neglecting terms, or by a combination of the three procedures. Thus, a solution of an 
approximation to the problem is obtained. 

Even though solutions may not be attainable to general transportation and production 
problems, a statement of a general problem can serve useful purposes. By comparing a 
simplified statement of a problem with the general statement, one can often determine the 
importance of omitted terms and the effect of approximating nonlinear functions by linear 
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functions. Also, examinations of a general statement may provide interesting qualitative 
results without an explicit solution being obtained. 

The following sections will include a statement of a general transportation and 
production problem, a discussion of the restrictions placed on the general statement to obtain 
a well-known problem, and derivations of several results from the general problem, even 
though the problem is not solved. 


THE GENERAL TRANSPORTATION AND PRODUCTION MODEL 
Transportation Rates 

Since in a continuing operation costs depend more directly on rates of production, rates 
of demand, and rates of transportation of an item than on quantity, these rates are taken as the 
fundamental variables. Let 


An; be the number of items of the k-th type transmitted per 
unit time from terminal n towards terminal j. 


It is assumed that Fin: is positive when goods are being sent out of termina] n along 
the transport path (n,j) and is negative when goods are being received at terminal n along 
that path. The algebraic flow Pin of the k-th type of item through terminal n is represented 
by 


Plan = nj 


If Pin is negative, then the rate of receiving exceeds the rate of transmitting goods of the k-th 
type at terminal n. If at terminal n there is a known demand rate Win >.0, then we may 
specify 


Pun = - Wen Win = ~ 4 ins 


If Pin is positive, then the rate of transmitting exceeds the rate of receiving goods of the k-th 
type at terminal n, and to keep the system in balance a supply rate of magnitude Pin must be 
maintained at that terminal. One way of meeting a supply rate, Pin ” 0, at terminal n is by 
having a production facility at the terminal. 

An unknown demand may occur at a terminal. For example, if at terminal n, goods of 
the k-th type are converted to goods of the (k+1)-th type, i.e., 


Pin = Pkn Px tn) ? 


where Paid n is a supply rate necessary to keep the system in balance, then Pin < 0 is the 
demand rate of the k-th type of goods necessary for the conversion process, which is to 
provide the supply rate of the (k+1)-th type of goods. The preceding relation may be writtenas 


. Fnj = Pkn th I+ 1.) , 
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where Ly represents the sum over transport paths terminating at terminal n for the k-th 
type of goods and the Lo represents the sum over transport paths terminating at terminal n 
for the (k+1)-th type of goods. 


Losses 
On each route between terminals the loss or breakage rate is considered to be a known 
function of the transport rate over that route. Let 


Linj be the number of items of the k-th type lost per unit time 
in transporting Inj items per unit time from terminal n 
towards terminal j. 


The rate at which items are sent over a route less the rate of loss must give the rate of 
receipts. Hence the loss equation is 


Any ~ Linj = ~ gn - 


The loss equation defines Finn? for 


Ginn = Linn’? » 


and — is the breakage or handling loss at terminal n. The loss equation may be used either 
as a set of further relations among the transport rates or as a means of replacing about half of 
these rates, say those Inj with n>. j, by expressions in the others. 


Costs 


The cost of operating a terminal is considered to be a function of the flow through the 
terminal. Let 


Fin be the part of the cost per unit time of operating terminal 
n allocated to the k-th type of product. 


The Fin may be considered as a function of Pin (or, equivalently, of Ln; or Win)» or it 
may be considered as a function of the Vnj for all j. 

The cost of operating a transmission line or route between terminals is taken to be a 
function of the traffic over that route. Let 


Tynj be the cost per unit time of moving Gin items of the k-th 
type per unit time from terminal n towards terminal j,. 


By defining Ten in this fashion, we may separate costs of sending and costs of receiving 
from each other. If there are no costs involved in the receiving of items, then 


Tiny = ° 


for nj <0. 





G. W. EVANS II 


The total cost to be minimized is 


Stockpiling 

A stockpile of a specified amount desired by a certain time at terminal n may be 
interpreted as a stockpile rate Sin The stockpile rate may be incorporated in the model in 
several ways, depending upon the manner in which one desires to account for its costs. It may 
be set up as a demand rate at a new terminal, n', connected to terminal n by a single trans- 
port path (n,n'). This arrangement would allow the assignment of separate loss rates Linn’ 
to the handling of the stockpile and Lin'n' to deterioration of goods held in the stockpile and 
would allow the assignment of separate cost rates Tue! A Tin'n? and Fin’ The stockpile rate 
may be considered as a part of the demand or production rate at terminal n. That is, if Win 
is the known demand rate at terminal n for the k-th type of goods and Sin is the desired 
stockpile rate, then we may write 


W,, = W * Sn = 72 Ung - 


Similarly, if Syn Occurs at a terminal not having a known demand rate, we may write at that 
terminal 


Phen = Sin + nj 


In this case, the cost associated with Sin is incorporated in Fin’ Another method of intro- 


ducing Sin would be by including it in the term Ann? 1es if the loss rate at terminal n is 
Linn’ then we may write 


Qnn = Un’? = xan + Sqq)/2 - 


In this case, the costs associated with Syn are incorporated in Tynn’ Thus, in the following 
discussion we will assume stockpiling to be accounted for in one of the preceding ways and not 
refer to it by special symbols. 


Restrictions 


Very often physical restrictions play an important role in the statement of the model. 
Restrictions on the rate of transport over a particular path, such as 


0< Max || dyn! ’ | icin |} < ny? 


which restricts the maximum rate of transport over any part of a path, or if the flow over a 
path is to have a certain direction, say from n to j, and to be less than a maximum rate, then 


inj + Linj < Akin < O< Vcnj S Qenj 
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may be imposed. Eliminating a transportation path for a given type of product is equivalent to 
setting Ovni = 0 and Lanj = 0. Finally, upper bounds, Puy to the production rates at sources 
are to be expected: 


O< Pun < Pin - 


If a solution is to exist to a problem having these restrictions, then there are obvious 
physical constraints on the system which must be satisfied, such as 


LZ Pin2 pi Win » 


Mp By 


for each type of goods, where ): represents a sum over terminals having factories and 


represents a sum over terminals at which demands occur. Where maximum rates of transport 
over a path are involved, then 


for each demand. In fact, by having restrictions on the amount of production or on the maxi- 
mum transport rates or both, demands can be increased until they cannot be met, and no 
solution will exist to the problem. For the problems discussed here, it is assumed that a 
solution exists. 


COMPARISON OF A LINEAR TRANSPORTATION MODEL WITH A GENERAL MODEL 

In this section, the Hitchcock-Koopmans [1,2] transportation problem is derived from 
the preceding general transportation and production model to demonstrate where terms may 
have been linearized or omitted. Evaluation of the effects of these simplifications must be 
made with respect to specific problems and is omitted here. 

A Hitchcock-Koopmans model is of a transportation problem in which a homogeneous 
product is shipped. In terms of the definitions of the preceding section, the subscript k 
assumes only the value 1 and is omitted as a subscript in this section. The allowable trans- 
portation paths are only those between origins and destinations. The origin terminals are 
designated by i<n<ny - 1 and the destination terminals by ngs j<N; the subscripts n and 
j may be considered as having only these ranges. No losses are considered to occur along the 
transportation paths. In fact, losses are not considered as a part of the model, i.e., 


day = “Mn 


Thus, the rate of shipping at one end of a transportation path is the rate of receiving at the 
other end. The shipping rates at the origins as well as the demand rates are specified, 
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for 1<n<n,-1 and 


an} an} 
-p, = Wy =-) | din = 2. Gyj 2° 
n=1 n=1 


for ny<j<N, where the Wh and W; are known. The restrictions on the model are that the 
movement of goods be only from origins to destinations, 


for l<n< Ng - 1 and mgs j<N and that the total of the shipping rates from the origins be 
balanced by the total of the demand rates at the destinations, 


ng -1 N 
dD. %,* ‘e ”" 
n=1 j=ny 


Costs for operating the terminals are omitted, 


F, = 0 


for 1<n<nj,-l, and 


tae 


for Ny < j < N, leaving the shipping costs as the only ones considered; these costs are assumed 
to be linear, 


Tay = ©nj Inj ’ 


where the Cay are known constants for each transportation path. 


The desired solution of the Hitchcock-Koopmans model is the set of Qj that satisfies 
the conditions discussed above and that minimizes the total cost of shipping: 


Ny -1 N ny -1 N 
c=) De Ty" 2, 2 cate 
nm=1 j=ny n=1 j=ny 


The solution of this problem may be obtained by linear programming methods such as the 
Simplex Method of G. Dantzig [3]. 
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USE OF TWO-PRODUCT SYSTEM TO OBTAIN RESULTS FOR A MULTIPRODUCT 
SYSTEM 

In this section, transport paths of the general transportation and production problem 
are restricted in such a manner that results concerning multiproduct systems may be derived 
from a two-product system. The number of transport paths are restricted in the following way. 
The terminals 1< n <N are divided into I-groups ny < n<n, 4-1, 1<i<I, where n= 1 and 
my, 471 = N. Only transportation paths (n,j), from terminal n to terminal j are allowed which 
satisfy 


my Shs my ge 


ns3<™,27) 


for 1< i< I-1. For example, if N= 6 and n, = 1, ny = 3, ng = 5, and n, -l=n,,,-1=6=N, 
the allowable paths are: 


(1,1) (2,2) (3,3) (4,4) ~««(5,5) (6, 6) 
(1,2) (2,3) (3,4) (4,5) ~—«5, 6) 
(1,3) (2,4) + (3,5) ~—«(4,6) 

(1,4) (3,6) 


The paths not allowed are (1,5), (1,6), (2,5), and (2,6). Under these restrictions, results which 
are extensible to larger systems may be obtained from systems involving three groups of 
terminals. An example of such a system involving three groups of terminals is the following 
two-product system. For this system, we will assume a minimum cost solution exists which is 
obtainable by the method of Lagrange multipliers. For a solution to exist by this method 
implies that the restrictions on maximum production and transport rates are automatically 
satisfied and may be omitted from the discussion, If, by the method of Lagrange multipliers, 
the minimum cost solution exists and satisfies the physical constraints, then it is the desired 
solution for the system. But should this minimum cost solution exist and not satisfy the 
physical constraints, it indicates where the system should be operating, i.e., what additional 
production and transportation facilities should be added to the system. 

In a two-product system, products of both types may be transported over the same 
path. For example, one might conceive of a paper industry in which paper stock is made from 
wood pulp at the terminals for which 1 < n < Ny, finished paper is manufactured from paper 
stock at the terminals for which n+ 1<n<nj-l, and demands for finished paper exist at the 
terminals for which ny <n<N, In such an industry it may be necessary to transport paper 
stock from a paper-stock factory to the desired finished-paper factory via one or more other 
paper -stock or finished-paper factories because of economical routings for the system. 
Similarly, finished paper transmitted to a demand may have to be routed through several 
terminals. A possible result of such movements is that both paper stock and finished paper 
would be transmitted along the same path. The model for the system may be constructed as 
follows, where lower-case and upper-case letters are used to distinguish the quantities 
associated with each product and the subscript k is dropped. 




























Qj 


nj 


nj 


system: 


(1) 


be the rate at which paper stock is produced (or demanded) at terminal n, 


1<n<njy-l; 


be the cost per unit time for operating terminal n assigned to the paper- 
stock phase of the system, 1 <n< ng-1; we assume f, = f,,®,); 


be the rate at which paper stock is transmitted from terminal n towards 
terminal j, 1 <n,j< ny-1; 


be the cost per unit time for transporting paper stock at the rate yj 
from terminal n towards terminal j, 1 < n,j < nj-1; we assume 


tj = taj nj)» n fj, and tan = 0; 


be the loss rate of paper stock in transporting paper stock at the rate qj 
from terminal n towards terminal j, 1 < n,j< nj-1; we assume 
lj = nj Inj)» jn, and 1, is known and constant for each n; 


be the cost per unit time for operating the paper-stock phase of the 
system; 


be the rate at which finished paper is produced (or demanded) at terminal 
n,n, + 1<n<N; 


be the cost per unit time for operating terminal n assigned to the 
finished-paper phase of the system, a, + 1 <n< N; we assume 

be the rate at which finished paper is transmitted from terminal n 
towards terminal j, ny + 1 <n,j <N; 


be the cost per unit time for transporting finished paper at the rate Qi 
from terminal n towards terminal i, Ry + 1<n,j <N; we assume 


Thi - Tyj(Qnj)» nf j, and Ta = 0; 


be the loss rate of finished paper in transporting finished paper at the 
rate Qnj from terminal n towards terminal j, ny + 1<n,j <N, we 
assume Lyj = Lj (Q,)) n # j, and L,,, is known and constant for each n; 


be the cost per unit time for operating the finished-paper phase of the 
system; and 


be the cost rer unit time for operating the complete system. 


In addition to the above definitions, the following restrictions are assumed for the 


P,=P,®,) OF Py=PAP,), m+ l<n<n-l, 
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i.e., the conversion rates of paper stock to finished paper are known; 


(2) Pi -Wn> nycn<N, 

where the Wi >. 0 are the known demand rates; 

(3) pj ~ Any = ~94n > l<n<n-l, n<j<n,-l, 
and 

(4) Qj 7 Lay = Qn n,+i<n<N, n<j<N, 


i.e., for each type of product there is one loss equation for each transportation path and 
terminal; 


ng-l 
(5) P=) Wy: 1<n<m); 
1 
ng-l 
(6) p,?,) = >. Any » n,+ l<n<ng-l; 
j=1 
N 
(7) P= >; Qyj> Mpt tsnsng-l; 
j=n,+ 1 
and 
N 
(8) -W, = >; Qj? mygcn<N. 
j=n,+ 1 
Finally, the total cost per unit time of the two-product system may be written as 
(9) y=zec+C, 
where 
ny-1 ny -1 ny -1 
(10) c=) t0,)+ DD. tan) 
n=1 n=1 j=l 
and 


N 


N N 
> Fat 2, 2, Tay@y)- 


=> 1 n=n,,+ 1 j-n,+ 1 


(11) ¢ 
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The desired solution is the set of q,. and Qn 
straints of Eqs. (3), (4), (5), (6), (7), and (8). 

Generally, for the two-product system, the a cannot be determined in terms of the 
Qa and one cannot assume a single type of product, for example, the dollar value of the 
product, to flow through the system. To show this, Eq. (1) is written as 


that minimizes Eq. (9), subject to the con- 


ng N 
(12) P Gy = Py,) = Pa im Qaj} > Bp * l<n<nj,-l. 
j=1 jen,+1 


There are (ng-n,,-1) equations of this type relating qyj'8 with Qj * Furthermore, some of 
these q,;'8 are related by (ng-n,,-1) [1+ (ng-n,,-2)/2] additional equations, Eq. (3) where 

n, + 1<n<nj-l and n<j<nj-l. Thus, there are (ng-n,-1) [2 + (ng-n, -2)/2] equations to 
determine (ng-n,,~1) (ng-1) quantities qj: For the qj to be written as functions of the Qj 


we must have 
(ng-n,-1) (ng-1) - (ng-n,-1) [2+ (ng-n,, -2)/2] = (ng-n,-1) [n,-3 - (ng-n,-2)/2] <0. 


Since by definition ng-l >, ny >. 1, the preceding expression is zero only for the smallest 
allowable values of n, and Ny» Ng = 3 and n= 1, i.e., when there is only one paper-stock 


factory and one finished-paper factory. The implication is that when there is only one paper 
stock factory and one finished-paper factory, connected by a single transport path, the two 
factories may be considered as a single factory making finished paper. Otherwise, the expres- 
sion is greater than zero and the qj cannot be determined in terms of the Qn 


The preceding result may be generalized to systems containing more than three groups 
of terminals, since any group of terminals may be considered as demands for the preceding 
group. For transport-path restrictions of this section, the general result is: any two groups 
of terminals, each consisting of a single terminal, which are connected by a single transport 
path, may be considered as a single terminal producing the product of the latter group. 
Heuristically, this is a trivial result. Not so obvious, however, is the result that other 
combinations of transport paths and terminals do not allow a reduction in the number of types 
of products. 

To obtain a stationary solution for the two-product system by the method of Lagrange 
multipliers, let the independent variables be dy 1<n,j < ny-l, Qa n+ 1<n,j<N, and 











Py ny +1<n <ng-l. Define 
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n=1 n=n,+ 1 n=l] 


n, ny-l ng-l ny-1 ny-1 
ale >. Ap «| ” ta fn [Pa Pa] . } at taj nj) 
n=1 os" 





ny-l N 


N N 
+ D° FA@)+ pie F,(-W,) + 2. ¥ T,y(Q,y) 
n=n_+1 n=n nm=n+1 j=n.+1 
n, d "pjyn P 
ng-l ng-l ny-l N 
7 z By Pay) - ) Qj} ~ ym M, | Pa - pM Qj 
n=n,+ 1 j=1 n=n,+ 1 j-n,+ 1 
N N Mg-2 Ny-l 
+ 27 M(Mat 2) ay) - 22 2) rng ng + Syn - Aap) 
n=n, jen,+1 n=1 j=n+1 
N-1 N 
(13) — Dy QE Sng @ny + Yn ~ Lng) 5 
n=n,+ 1 j=n+1 


and set 0Z/8 d,, = 0, 1<n,j <ny-1, jfn; 02/2 Q,; = 0, ny + 1<n,j<N, jfn; and aZ/0P, =0, 


S- | n+ l<sn<ny-l. That is, for 1<m<n, and n+ 1< j<nj,-l, 
1S (14) az Mn yy Sy = 
°4nj dp, dq, j nl dq,j 
For 2<n<m, and 1<j<n-l, 
df. dt. 
, (15) 02 nL n= 0- 
24nj Fy My 


For n, + 1<n<nj-2 and n+ 1<j<n,-l, 


ny py 
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dt 
(16) oo ana “Zino, 
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For n, + 1<n<nj-l and 1<j<n-l, 


az Sy 
2dnj Mn 





(17) 





+ Hy ~Agn = 9- 


For n+ 1<n<N-l and n+1<j<N, 








aT dL 
(18) =. ah. M,- Ay -_ =0. 
Qj FRaj j 
For n,+ 2<n<N and n+ 1<j<n-l, 
az ST yy 
(19) 30 AQ eB On 93 


and for n+ i<n<n,-l, 


df_\ /dp.\ dF dp 
(20) 02 _(_2 —*)\+—*-y(— |-m,=0. 
oP, \dp,/\aP,/ dP, dP, 


Equations (3) through (8) and (14) through (20) are used to determine the values of Gy? Quay 
Any Any Hy and M, (note q,, =1,,/2 and Q., = L,,,/2 are determined from Eqs. (3) and 
(4), respectively, when n= j), which represent a solution of the two-product system. 


The solution of the two-product system may be considered as the solution to three 
related problems. Suppose, for example, the finished-paper phase of the system is considered 
separately. If a solution is desired as a function of the parameters P_, ny + 1<n<n,-l, 


then the solution for this subsystem is obtained by defining a quantity ¢ as 


nq-l N N N 
g=c- )° m,(P,- ). Qy)* 2, M, (W, + 7 Qnj 
n=n,+ 1 }-n,,+ 1 n=n,y j=n,,+ 1 
N-1 N 
— Di Qa Sag Qag + %m - Uy 
n=n,,+ 1 j=n+ 1 


and setting the 0¢/0Q,, = 0 for n, + 1<n,j<N and jfn. This implies that Eqs. (18) and 


(19) must be satisfied. Functionally, then, Eqs. (4), (7), (8), (18), and (19) may be solved for 
the Qay M. and Aaj as functions of the P,, and known constants; i.e., 
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Q,j = Qj (Poy errs Png-t) 


(21) M, = M, (Pays yosoes Pat) 


Any =Anj (Paget Png-t) ; 


Similarly, if the paper-stock phase of the system is considered alone and if a solution 
is desired as a function of the Py the solution for this subsystem is obtained by defining a 
quantity e as 


ny-1 nq-l ny-2 ny-l 
nines Z Hy| PyP,) - Fa yj | - th R nj nj + Fyn ~ tng) 
n=n,+ 1 j=l n=1 j=n+l 


and setting the de/ddy; = 0 for 1 <n,j < ng-1 and j#n. Thus, Eqs. (14), (15), (16), and (17) 


must be satisfied. Functionally, Eqs. (3), (5), (6), (14), (15), (16), and (17) may be solved for 
the qj? My» and raj as functions of the Py i.e., 


Qj = Qj (Pa,+1 eee 
(22) My = Up (Page aus Paarl) 


aj = nj (Pa, geeees Pag-t) 


Finally, the finished-paper factories may be considered as a subsystem. The solution 
is obtained by defining the quantity E as 


ng-l ny-l ny-l N 
=e >. F,(P,) + b i fn [Pan ~ Z. M, P,- Pa Qn 
n=n,+ 1 n=n,+ 1 n=n,,+ 1 j-n,+ 1 
ny-1l ny-1 
7 Zz Hy | Py(P,) - z. Qj 
n=n,+1 j=1 


and setting the dE/oP, = 0, n, + 1 <n<nj,-1l. Thus, Eq. (20) must be satisfied. Functionally, 
Eq. (20), after the dy? Qa’ Has M,, Any? and ay have been replaced by Eqs. (21) and (22), 
may be solved for the PY n+ l<n< ny-l. 


Thus, the two-product system may be made to appear as three simpler problems. That 
is, the same equations are involved whether the two-product system is considered as a whole 
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or is divided into three parts. For a system containing I-groups of terminals having the 
transport path restrictions of this section, the preceding result becomes: an (I-1)-product 
system may be considered as I-1 single-product transportation and production problems 
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related by I-2 production problems in each of which one product is converted into another. 





ON OBTAINING A SOLUTION TO THE TWO-PRODUCT SYSTEM BY THE METHOD OF 
LAGRANGE MULTIPLIERS 
As in the preceding section, we assume that a solution exists for the two-product 
system and may be obtained by the method of Lagrange multipliers. Thus, Eqs. (3) through (8) 
and (14) through (20) represent a valid system of equations to be solved for the unknowns dy? 


Qy Any A nj’ Has and M,- In this section we will show that the number of equations may be 


greatly reduced with a corresponding reduction in the number of unknowns. 
The range of values for n and j in Eq. (19) is n, + 2<n<WN and n, + 1<j<n-l, 


which may be written as n+ 1<j<N-1 and j+ 1<n<N. By interchanging the roles of n 
and j, Eq. (19) may be written as 


(19a) 


where n, + 1<n<N-l and n+ 1<j<N. Substituting the right side of Eq. (19a) for Anj in 
Eq. (18) gives 


(18a) 


and n, + 1<j<n,-l and by Eq. (14-2) for the range 1<n<n,-l and n+ 1<j<n), In 
Eq. (15) the range of values for n and j is 2<n<n, and 1< j <n-1 which may be written 
as 1<j<n)-l and j+ 1<n<n, By interchanging the roles of n and j, Eq. (15) may be 


written as 


(15a) 


where l<n <n,-l and n+1<j <M). Substituting the right side of Eq. (15a) for d nj in Eq. 
(14-2) gives 





is 


for n, + 1<n<N-l and n+ 1<j<N. Thus, the unknowns ‘yj have been eliminated, and the 


two sets of equations, Eqs. (18) and (19), have been combined in Eq. (18a). 
The range of values for n and j in Eq. (14) is l<n<n, and n+ 1<j<nj-l, which 


may be divided into two ranges. Let us designate Eq. (14) by Eq. (14-1) for the range l<n<n, 
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df, dt, /df 
(14-2a) ee, ee 
dp, dqj dp; ddin dq,j 


for 1<n<n,-l and n+ 1<j <n). 

The range of values for Eq. (17) is n, + 1 <n<ng-l and 1<j<n-1, which may be 
divided into two ranges. Let us designate Eq. (17) by Eq. (17-1) for the range n, +1<n< ng-l 
and 1<j<n, and by Eq. (17-2) for the range n, + 2<n< ng-l and n+ 1<j<n-l. By 
interchanging the roles of n and j, Eq. (17-1) may be written as 


(17-1a) r tin 
-la 24+ S 
nj 43 * aq. 
where 1<n< ny, and n, + 1<j <ng-l. Substituting the right side of Eq. (17-1a) for nj in 


Eq. (14-1) gives 
df dt dt dl 
(14-1a) ame + _w - My + _in _ =0 
Pn “nj jn) \ ny 
for 1<n<nm, and n+ 1<j <Ng-l. The range of values for n and j for Eq. (17-2) may be 


written as Ny + 1<j <Ngq-2 and j+1<n <ng-l. By interchanging the roles of n and j, 


Eq. (17-2) may be written as Eq. (17-1a), which will be designated as Eq. (17-2a) for the range 
ny + 1<n <ng-2 and n+ 1<j<nj,-1l. Substituting the right side of Eq. (17-2a) for nj in 


Eq. (16) gives 


di dt dl 
(16a) be My - uj+—= _ a =0 
aq, j dain dq,j 
for n, + l<n <Mg-2 and n+ 1<j <ng-l. Thus, the unknowns Anj have been eliminated, and 


the sets of equations, Eqs. (14), (15), (16), and (17), have been combined to form Eqs. (14-1a), 
(14-2a), and (16a). 

Next, we note that the conversion rates of paper stock to finished paper are known and 
dp, /aP,, #0, even for P= 0. Thus, the unknowns H, may be eliminated by writing Eq. (20) as 


df dF dp 
(20a) mma (at) Kes 
n 


dP, 


where n, + 1<n<n,-1, and substituting the right side of Eq. (20a) for HM, and My in Eqs. 
(14-1a) and (16a). 
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The M, may be eliminated by using a subset of the equations of Eq. (18a). For 
example, we might choose n= N-1 and j = N and write Eq. (18a) as 


d 
dby-1,N_ ce . : | Ty N-1  9TN-1,N 


— ? 
dQy-1.N dQy-1,n/ Wy w-1 4@n-1,N 





(23) My-1 +( 


for n= N-2 and j=N, Eq. (18a) becomes 


. a aT aT 
(24) My 2 + eae _ ) My = . i ea) N,N-2 N-2,N E 


> 
dQy-2.N dQy-2n/ Ty n-2 I@y-2 Nn 





and for n= N-2 and j= N-1, Eq. (18a) becomes 


dly_2 N-1 ve dLy_2n-1\?Ty-1,n-2 9TN-2,N-1 
Go) Mes*lag i) MI a dQ dQ 
N-2,N-1 N-2,N-1/ 9@n-1,n-2 9@y-2,n-1 


Eqs. (23), (24), and (25) may be considered as three linear equations in the unknowns My_»2: 
My_;, and My and solved for these quantities since, in general, the 1 - dL/dQ #0. By using 
j=N and n, + 1 <n <N-3, Eq. (18a) may be written 


alin qT yn aT LN 
(26) M, = f ae) . + | = dQ. ; 


and the remaining M,> My + 1<n <N-3, may be obtained. Having obtained expressions for 
the M,, a, + 1<n <N, these expressions may be substituted for the M, and M; in Eqs. 
(14-1la) and (16a) and the remaining equations (those for the range n, + 1<n<N-3 and 
n+ 1<j<N-1) of Eq. (18a). Thus, the unknowns M,, may be eliminated and the number of 
equations of type Eq. (18a) reduced by N-n, equations. 

Finally, the substitution of PY as defined by Eq. (7), into Eq. (6) gives 








N ng-l 
(6a) Ph pa Qyj = Pe qj 
j-n,+1 j=1 


for n+ 1<n<n,-l. 

The basic system of equations to be solved are Eq. (14-la) for 1<n <0) and 
n,+ 1 <j <nj,-1, Eq. (14-2a) for 1<m<n)-l and n+ 1<j <M, Eq. (16a) for n,+ i<n<n,-2 
and n+ 1<j<nj,-1, Eq. (18a) for n, + 1<n<N-3 and n+ 1<j< N-1, Eq. (6a) for 
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ny + 1<n<n,-l, and Eq. (8) for ny <n<N. It is assumed that, in this system of equations, 


the 1, and M, have been eliminated as described above, the p,, have been eliminated by 

Eq. (5), the P, have been eliminated by Eq. (7), the q,; for n >j have been eliminated by 

Eq. (3), and the Qj for n >j have been eliminated by Eq. (4). Thus, the basic system of 
equations are solved for the nj l<n< Nq-2 and n+1<j< ny-l, and the Qj nyt 1<n<N-l 
and n+ 1 <j <N, which is the minimum number of unknowns in the two-product system. 


CONCLUSIONS 

Since this paper deals primarily with the statement of a problem rather than its solution, 
it is appropriate to include some remarks concerning this form of statement and a discussion of 
the types of solution desired. 

With respect to the statement of the problem, we consider the rates rather than the 
quantities of goods being produced, transported, or converted. This choice was made because 
orders for goods may be best considered in terms of demand rates. That is, the promised 
delivery date of a given quantity of goods prescribes the minimum rate at which those goods 
must be produced and transported, after lag times have been deducted, in order that this 
demand may be met; and, thus, the demand should be represented as a rate so that its effect 
upon shipment losses will be reflected in the solution. 

If a stockpile of a given type of goods is desired or is to be depleted at a given terminal 
by a given date, then a demand or production rate, respectively, is specified at that terminal. 
Also, the rate statement of the problem allows steady-state solutions without carrying time as 
an independent variable in its statement or solution. 

Goods introduced into the system to help meet a demand Win where the cost of 
production and of transporting these goods are not a part of the system also may be included 
in the general model. For example, if these goods are introduced along the transport path 


(j,n) at the rate Gy; - and along this path the loss rate is Lig then we may write 


Lin = Lin ~ Gain > 


and the loss equation for that path remains 


Gin ~ “ign = ~Ixnj - 


In this paper the method of solution discussed has been that of Lagrange multipliers, and 
constraints, such as maximum production rates and maximum transport rates were not considered. 
If these constraints are incorporated, then a direct-gradient method of solution, such as the 
one of Ablow and Brigham [4], may be employed. Such a method is used to solve a system of 
constraints in the form of inequalities, and each equation used in the general model to describe 
a constraint (such as a loss equation) would be replaced by two inequalities. This can be done 
since f (x,, Xp, eves x) = 0 implies that both inequalities f (x,, Kg, +-+5 x) >0 and 


f (x4, Xo, +++, x) <0 must be satisfied. The solution, however, when it exists, by the method 


of Lagrange multipliers provides valuable information, even when that solution does not fall 
within the complete constraints of the system. It describes where the system should be 
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operated, and from this information one may determine where additional production and trans- 
portation facilities should be added to the present system. 

The preceding remarks hint as to the type of solution desired. Suppose the criterion 
function is minimized for a given set of constraints, and then the constraints are changed, that 
is, new demands are incurred or old demands are fulfilled, transportation paths are added or 
removed, or a production plant is put into operation or shut down. Rather than solving a 
complete problem, one would prefer to ask how the minimization of the criterion function is 
affected by these changes. In fact, this suggests that the method of solution desired should 
prescribe the change in the allocations of production that produces a new minimum of the 
criterion function caused by a single change of constraint. When there is more than one change, 
these changes should be programmed one at a time. 
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PROBABILISTIC ERRORS IN THE LEONTIEF SYSTEM* 


R. E. Quandtt 


Princeton University 


In order to judge the reliability of input-output estimates one has to consider the effect 
of errors in input coefficients on the solution of the Leontief system. The existing treatments 
of this problem customarily establish bounds for the solution or for the elements of the inverse 
matrix on the basis of various assumptions about the magnitude and the direction of the errors 
[1-7]. With minor exceptions the previous treatments of the problem have not considered the 
possibility that errors in input coefficients may have stochastic properties [8]. This article 
explores some of the implications of the assumption that input-coefficients are distributed 
according to some probability distribution. 


1. THE RATIONALE FOR THE ASSUMPTION OF PROBABILISTIC ERRORS 
The Structural Rationale 

If it is assumed that demand curves fluctuate in a probabilistic manner [9], i.e., that the 
quantity demanded is a random variable, it is relatively easy to demonstrate that the factor 
proportions employed by the firm and hence the observed input coefficients themselves are 
subject to random changes, provided that some of the assumptions of input-output analysis are 
relaxed. Observed input coefficients will tend to be distributed according to some probability 
distribution if at least one of the following statements is true: 

(a) Firms producing identical commodities have identical production functions which 
are linearly homogeneous and do not allow factor substitution, but at least two distinct activi- 
ties are available for the production of a given commodity. 

(b) All firms producing identical commodities have identical production functions, and 
factors are strictly complementary, but the firms' expansion lines are not straight lines 
through the origin. 

(c) A given firm employs only one activity in order to produce a given commodity, and 
the production functions of all firms are linearly homogeneous, but different firms either have 
different production functions or face different relative factor prices. * 





*Manuscript received February 24, 1958. 

tlam indebted to Professors Wassily W. Leontief and William J. Baumol, Mr. Michio Hatanaka, 
and Dr. Hartley Rogers for valuable advice and criticism. Ofcourse, Ialone am responsible for 
errors. 

‘Note that these alternative assumptions are as weak as possible. In other words, we have re- 
tained in each of the alternative assumptions as many elements as possible that are consistent 
with the standard interpretation of input-output analysis. 
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If any one of these conditions is fulfilled, input coefficients observed at different times 
will generally be different. Furthermore, it is not meaningful to speak of "true" coefficients, 
because they would exhibit variations even if they were obtained by taking exhaustive samples 
every time. 








The Sampling Rationale 

Usable input-output tables are not very numerous, and the reason for their scarcity is 
mainly the mechanical difficulty of compiling them. The labor of collecting data and the costs 
of collection might be reduced by resorting directly to random sampling of coefficients.* One 
could (at least theoretically) obtain the value of each input used by a firm per dollar's worth of 
output by taking a random sample of individual firms stratified according to the desired indus- 
try classification. The sample mean for each coefficient would provide an estimate of the 
derived coefficient. This method would not obviate the necessity of giving careful attention to 
the usual types of conceptual and computational problems arising in the construction of input- 
output tables: the distinction between producer's and purchaser's values; the treatment of 
special categories such as inventories, gross capital formation, government, foreign trade, 
households, etc. The unallocated sector would theoretically disappear, and its actual existence 
could result only from an individual firm's inability to allocate its expenses to nondummy cate- 
gories. The crucial difference would be that one would not generally be able to derive a flow 
chart consistent with independently obtained control totals. 

If such an approach were employed, one would be interested in the following questions: 
(a) What is the meaning of the solution of such a probabilistic system? (b) Can one theoreti- 
cally attach confidence limits to the solution? (c) Is it possible to calculate the moments of 
the distribution of the solution? (d) How are the moments of the distribution of the input coef- 
ficients related to the moments of the distributions of either the elements of the inverse or of 
the solution? (e) What, if anything, can be said about the distribution of the input coefficients ?* 

The basic tool of input-output analysis is a system of equations, 


(1.1) ee aij * = Yi Tewrrr. 3 
j=l 


where the x's 2 re the gross outputs in constant dollars of n commodities, the a,i'8 are the 
input coefficie its (defined as the fraction of a dollar's worth of commodity i needed as current 
input to prodr.ce one dollar's worth of commodity j), and the y's are the quantities in dollar 
terms of the various commodities not used up in producing commodities. The y's therefore 
represent "shipments to final demand" or the "bill of goods." In matrix notation, (1.1) is 
written as 





*Present-day input coefficients are obtained by certain indirect estimating procedures which 
in most instances are equivalent to a complete count. See W. W. Leontief, 'Some Basic Prob- 
lems of Empirical Input-Output Analysis," in Input-Output Analysis: An Appraisal, National 
Bureau of Econ. Research (Princeton Univ. Press, 1955), pp. 9-22 (especially pp. 10-14) and 
"Comment" by A. Henderson, pp. 22-23. 
tNote that the mean of the population of ajj's (for fixed i and j)—which we may call the "true" 
coefficient—will generally differ from that obtained by an exhaustive sample of interindustrial 
flows. Any existing input-output table contains a certain amount of aggregation, at least to the 
extent that firms within an industry are aggregated. The arithmetic mean of coefficients in a 
sample is not an unbiased estimate of a coefficient as traditionally defined. 

*I hope to treat the last two problems in greater detail in future publications. 
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(1.2) (I - A)X ” Y, 





where I - A is the matrix obtained by subtracting the coefficient matrix A from the identity 
matrix, X is the column vector of outputs, and Y is the column vector representing the bill of 


goods. 


It is a customary restriction to assume that the column sums of the a,;'8 are all less 
than unity. If the elements of the coefficient matrix A are obtained by some random sampling 
technique, there is no guarantee that all column sums will in fact be smaller than unity. But 
without defending this position too vigorously, it will be assumed in the following that the re- 
striction on the column sums is always satisfied. This implies, of course, that the distributions 
of the input coefficients are defined for the interval (0, k), where k < 1. 

We now consider the entire coefficient matrix, and therefore the entire Leontief matrix 
I- A, as a random drawing from a population of such matrices. Ann xn matrix may be con- 
sidered as a point in n*-dimensional Euclidean space, and the probability distribution is defined 
on the set of n’-tuplets satisfying the a priori restrictions on such matrices. It will also be 
assumed throughout that the bill-of-goods vector is given and known with probability one. The 
solution of the Leontief system is given by the vector X, which will itself be a random variable. 

In general, one cannot assume that the elements of the Leontief matrix are independently 
distributed. Errors are frequently compensating, at least so long as coefficients are obtained 
by taking exhaustive samples [10]. When coefficients are obtained by random sampling, they 
may be independently distributed. The restriction rai <1 for all j does not necessarily imply 


* 


dependence. The distributions of coefficients in a particular column may be such that even if 
all coefficients assume their extreme values, their sum will still be less than unity. For sim- 
plicity's sake, it will be assumed throughout that coefficients are independently distributed. t 


2. THE EXPECTATION OF MATRICES 
The Theorems below are obvious and are stated without proof. 


THEOREM 1: Givena matrix P with independently distributed elements, the 
expectation of P, [E(P)], is the matrix of the expectations of the individual elements. 


THEOREM 2: Let P and Q be two matrices for which P + Q is defined. Then 
E(P + Q) = E(P) + E(Q). 


THEOREM 3: If C isa constant matrix, E(C) =C. 





*We define a probability measure on the space of n x n matrices such that singular matrices 
and nonsingular, non-Leontief matrices form a set of measure zero. Since the transformation 
which carries a matrix into its inverse is continuous and one-to-one, there will exist a prob- 
ability measure in the transform space also. The set of Leontief inverses, and hence the set 
of solutions, will have probability distributions defined on them. (Note that the restriction that 
all column sums of the A matrix add up to less than unity guarantees the nonsingularity of the 
Leontief matrix. See F. V. Waugh, "Inversion of the Leontief Matrix by Power Series," 
Econometrica 18:142-154 (1950). 
{Since the methods to be discussed involve nume rous approximations, it could probably be as - 
sumed without too much additional distortion that coefficients are not independently distributed, 
so long as the degree of dependence between them is relatively small. If x and y are two ran- 
dom variables (two input coefficients), one would then wish to stipulate that |E(xy) - 
E(x) E(y)|<e, where e is a small positive number and E the exrectation operator. 
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THEOREM 4: If k is a constant scalar and P a matrix, E(kP) = kE(P). 






THEOREM 5: Let P and Q be two matrices for which the product PQ is defined. 
If P and Q are independently distributed, E(PQ) = E(P) E(Q). 


COROLLARY: Since P and p-l are clearly not independently distributed, 
we have generally (although not necessarily) 


E(P) E(P™!) 4 E(PP~}) = E(I) = I. 


THEOREM 6: Let T bea "true" matrix and TX = Y. Let the observed system be 
given by PX = Y, suchthat E(P) = T. Then the expectation of the solution of the random 
system PX = Y is not generally equal to the solution of the true system. 


PROOF: The true solution is xX, = rly, and the solution of the observed system 
is X = pl Y. In order to prove the theorem, assume the contrary. Then 


rl 


X, - E(X) y - E(P7)y 


(2.1) 


(rt! - E(p-)jy = 0. 


If (2.1) is to hold for all column vectors Y, the matrix in square brackets must be the null 
matrix. Then 


gl . g(p7}) = 0. 


It follows that 
(2.2) I = E(P) E(P~4). 


By the corollary to Theorem 5, (2.2) is not generally true. This proves the theorem. Extend- 
ing this to the Leontief system, we can state that the expectation of the solution of a random 
Leontief system is not necessarily equal to the solution of the true system. 


3. APPROXIMATIONS TO THE VARIANCES AND COVARIANCES 
OF THE INVERSE ELEMENTS 

In order for confidence intervals for the solution to be established, the standard devia- 
tions of the elements of the solution have to be known. In order to guarantee the existance of a 
solution for any random drawing of a Leontief matrix to be guaranteed, it is necessary that the 
elements of the error matrix H be distributed such that the determinant |I- A- H|#0. This 
assumption is made throughout the rest of this paper. 

The solution of the Leontief system can be stated as 


n 
m= 2 by, iei,...,M, 
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where the bj;'8 are the elements of (I- ay. Then, 





n n on 
2 
(3.2) var (x;) = »: var (b..)y." + > covar (b.., b.,)y.y 
- vod et jz uP P97 
k#j 


for all i. The principal problem in obtaining the variances of the x's is to obtain the variances 
of the elements of the inverse and the covariances of the elements of the inverse with all other 
elements in the same row. Since coefficients are to be obtained by random sampling, it will be 
assumed that the means and the variances of the input coefficients are known or can be closely 
Ye approximated. It is required to approximate the values of the variances and covariances of the 
1 inverse elements on the basis of the properties of the distributions of the original coefficients. * 
Every element of the inverse is a fraction in which the numerator is a cofactor and 
the denominator the determinant. In order to evaluate the moments of the elements of the 
inverse, one has to evaluate the moments of fractions. The formulas to be developed below are 
quite general in that they do not depend on the density functions of the variables in question. t 
It is presumed in the analysis that the moments of the cofactors and of the determinant are 
already known.’ The results of this section are approximations, although the approximations 
can be carried to any arbitrary order. The necessary approximation formulas will first be 
derived in general and then, for illustrative purposes, applied to a hypothetical 2 x 2 Leontief 
system. 


em 


Variance of a Fraction 
The remainder of this section draws heavily on the work of several authors [11-13]. 
Let x/y be a fraction of two random variables with arbitrary densities. Then 


on wo wal) - (age 


y 


1d- Denote the respective deviations from the mean by 


x - E (x) = €5 
(3.4) 
y - E(y)=¥. 


Then (3.3) can also be written as 
via- 


2 72 
- (3.5) oa (2) _ p(B) + E E(x) + 4 7 
is 7 (E(y) + w)* E(y) + v 


The first term on the right-hand side of (3.5) is 





*W. Duane Evans [1] has suggested one possible approach. An alternative is suggested here. 
‘Certain conditions, however, must be fulfilled for their applicability,as we shall see presently. 
*There is no theoretical difficulty in calculating these moments, although in practice the cal- 
culations will be quite cumbersome. Simplifications are introduced by the assumption of 
independence among the coefficients. Clearly, the larger the dimensions of the system, the 
more cumbersome is the calculation of the relevant moments. 
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‘ [E(x)? [1+ ¢/E@)* an ay ee 1%. sy” a 
[EW f1+v/Ey) E(y) E(x) EY) (fey? 


E(x)\? ay 37 em | sey? 
*\EO)) 1" FO" ew? ¢ (F0 E@EW) E(x) (EW) 4 

















(3.6) +( Cau. ¢ ie ---) 
[EWP (EWP EGY) [EWPIEW!? 
provided that v <1. Whether or not this condition is fulfilled cannot be stated a priori. 








E 
Since the determinant of any Leontief matrix is between zero and unity, the condition will 
always hold if the expectation of the determinant is greater than 0.5.* 
The second term on the right hand side of (3.5) is 


2 2 ‘ 2 2 
em - (2) e(1+ =<) (- 4 v - +) 
E(y)+ ¥ E(y) E(x) E(y) (Eq)? 


g(x) \? |. vv. ow 4 
-(E@\ |, (1 - 5% on) 4 (sia a ae 
( E(y) E(y) [E (y)] E(x) E(x) E(y) g(x) [E (v7 




















- 2 
(3.7) = (2) 1 - E(¥) + E() ery Se) _ _5¢9) . | ae 
E(y) E(y) [E (y)/? E(x) E(x) E(y) 


Neglecting all terms of degree higher than second in (€W), one obtains the following second-order 
approximations to (3.6) and (3.7): 








(3.8) p(E@+e? . any ie sE(v") 4 E(Ey) ; wa } 


since E(¢) = E(V) = 0. Similarly, 








(3.9) [ 5G) + el j (ze (1 _ 2B)  2E(eW ). 


E(y)+y E(y) (Ey) =) EW) 


Therefore, 





*The variable y is the determinant and is restricted to the open interval (0, 1). Therefore the 
value |y| = ly - E(y) | can never exceed E(y) if E(y) > 0.5. 
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(3.10) 


y E(y) P 


2 2 
var(*) = fed jaw , Se. s5eW) 
(EY) [E@}? FE EG) 


and, by use of (3.4), 


tt) var (2) = (EO) (BOD, BOO) _ 2Bter) _) 





The Covariance 
We require the covariances of elements of the inverse with other elements of the 
inverse in the same row, i.e., the covariances of random variables x/y and z/y. Define 


a ome (G5) "5 -F) 05) 9G). 


Denote the deviations from the mean by 


x- E(x)=€, 
(3.13) y- E(y)=¥y, 
z- E(jz)=2, 


¥ ) and consider the first term on the right-hand side of (3.12). We have 








B(S.2)e8 (E(@) + EXE(2)+ 0) _ pone (+ ¢/E() (1+ el 
yy (EQ) +¥) [Eq (1+ y/E()? 





, 5. a ee ee eae 
“= ” Ewe E(@) * E@)* E@)E@)/\" EO) * eg? 














= EOEO | (1+ 5+ tye an Oe A 2 22 
x [E(y)}7 E(x) E(z) E(x) E(z) E(y) E(y)E(x) Ey)E(z) E(y)E(x)E(z) 
3y2 sy%¢ sy"z sy7ez. 
(3.14) +( 3+ 3 + 3 + )-- 
[E(y)}? (EQ) E@) (EQ)? E(2) [E()]}* E(x) E(2) 


provided that = s <1. Remembering that E(é) = E(¢) = E(¥) = 0 and neglecting all terms of 


E(y) 
degree higher than second in (éZ¥), one obtains 





2 
_ x 2) _ E(x) E(z) E(ét) — 2E(¢W — 2E(cW) | 3E(v') | 
O48) ETT) et [2 EG) EG EG) EG) E@ EO)” fegy? 


he 
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One can easily obtain from (3.7) the approximations 











x)_E@)[,, EW*) E(w | 
~_ *(5) "20 | * wt FRED 
Zz E(z) , E(w?) E(t y) 7 
3.17 ms a ; 
— (5) " #0) | * req?” EO EO | 








Toa second order approximation, 


x) _/z\ E(x) E(z){ 2EW*) EEW E(w) 
gay (5) B(5) - EWE f Tew?  E@EO SOV EGT | 








Then the covariance can be written as 




















covar (2,#)= E(x) E(z)[ Ev”) | Ee) BE EEW) 

yy) ey) [fem EO E@) EE) Elz) EW) 
or as 

2.19) covar(2,2) _E@E@)[ EQ") | Ez) EGY) Ezy) 

: y’y/ [E(y)]? [EWP E(x) E(z) E(x) E(y) E(z) E(y) | © 





4. APPLICATION TO THE LEONTIEF SYSTEM 

In order to make the calculations easier, the results of the previous section will be 
applied to a hypothetical 2 2 Leontief system. It will be assumed that (a) errors in the coef- 
ficients are distributed independently of each other, (b) the errors have zero means, (c) the 
density functions of the errors are symmetrical. The observed matrix is I - A - H, where H 
is the matrix of errors, or 


1- a), - hy, - 842 - hyp 


~ 89, - hoy 1 - ago - hoo 


(4.1) 


The assumptions state that E (h;;) = 0 for all i and j and that E(hj; h,,) = 0 for all i# r and 
j 7s. We shall also use the following notation: 


Diy Die 
(4.2) a- a)! = , 

boy boo 
(4.3) j1-A| =D, 


(4.4) 





|1-A-H| = 















(4 


fo! 
tel 


(4. 


ev 


(4. 


(4. 


wh 
the 


(4. 


(4. 


(4. 


(4. 


(4. 






ef- 








PROBABILISTIC ERRORS IN THE LEONTIEF SYSTEM 


for all i and j, where 6;, = 1 when i= j and zero otherwise. Finally, by rewriting (3.16) in 
terms of the definitions of the deviations from the mean, we obtain 





2 
(4.6) E(>) = E(x) |, , EW’) __E(y) 


y)  E(y) [Ey)}? EF) EW) 
The Expectation of the Inverse 
In order to obtain an approximation to the expectation of the inverse, (4.6) has to be 


evaluated for every element. First, 


(4.7) E(D*) = 


E(D**) = D? [1+ bz, E(h2,) + bey (hey) + be, E(hzo) + bt, E(h5,)] 





(4.8) 
+ E(h2,) E(hgy) + E(h25) E(hg,) , 
42 .. of 
(4.9) E(D*) _D'+DK+R_,,,,R. 
[E(D*)? D 


where K stands for the square-bracketed term on the right-hand side of (4.8) and R stands for 
the sum of the last two terms on the right-hand side of (4.8). Furthermore, 


(4.10) E(cg9 D*) = D[b,, D+ bop E(h3,)] , 
(4.11) E(c,, D*) = D[by9D+ by, E(n?,)], 
(4.12) E(-cy9 D*) = D[b,)D - by, E(hzy)] , 
(4.13) E(-cg, D*) = D[by, D - byp E(hg,)] - 


By the substitution of (4.7) to (4.13) into (4.6), 


— 


b. 
R 22 
by 1+K+—~ a “By Bh) Pia |) 1+K+—~ = by bb Eta) 


(4.14) E(I-A-H)7!= 
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The Variances of the Inverse Elements 
The application of (3.11) gives the requirea variances in analogous fashion. The only 
term which is yet to be calculated is E(x’). Calculating this for each element of the inverse, 


we obtain 






(4.15) 


(4.16) 


(4.17) 


(4.18) 
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E(c2,) = b#,D* + E(n2,), 


var (b,,) = bi 


var (by9) = bee 


2 


var (by,) = oa 


The Covariances by Rows 


E(c3,) = be D + E(hs,) ‘ 


Substituting the necessary terms in (3.11), we obtain 

















a ‘ . 
R  Elhgg) 2bo2 
K + — + 2 - E(ho9) 
D2 be, pe 4D : 
= 
a R E(ht,) 2b, 2 i 
i p? pe, vp? ae ; 
i R E(hi9) bo, og a 
K + — + pa a D2) 
c p* ow s é 
i E(h5,) 2big sg i 
hae be ee E(h9;) 
De ao n” 





square brackets in (3.19). We have for the first row of the inverse 


2 
E(xz) 941 942 DP 





E(x) E(z) _ bi Pas pe 


and the same for the second row. Then 


b 
R 22 
(4.19) covar(b,, bio) = Diy bie K+ — - 


pe 4D 






2 
E(hg9) + 


=1 


boy 
bj9D 











~~ 


In order to find the required covariances, we must still calculate the second term in 


E (bY) 
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(4.20) covar (bp bo) = b bp eee Pa E(hy,) + —? ad . 
. 1°22 21 "22 « e - i by D 1 


5. CONFIDENCE INTERVALS 

Various alternatives are open. A sample Leontief matrix could be obtained by ran- 
domly choosing a value for each of the n? aj;'8 from the population of coefficients. If we take 
repeated samples of Leontief matrices, the mean solution and the standard deviation of the 
mean solution can be calculated. Confidence limits for the mean solution are given by 


- ks < m< X¥ + ks, 


where X is the sample mean, s the standard deviation of the mean, m the population mean, and 
k a parameter derived from the t-distribution such that 


Pr[|X - m| > ks] <a, 


where ais the desired level of significance [14]. 
If the population variance is known, confidence limits for the solution (rather than the 
mean solution) can be found by Chebyshev's inequality: 


oA 


Pr[|x - w| >t] <—, 
t2- 


where x is any random variable, 1 its mean, and o? its variance [15]. This will generally 
result in wide estimates and is not recommended as an efficient technique. 

A confidence region for the solution (rather than for the mean solution) can be found 
more directly in the following manner. Let the density function of the coefficients of the input- 
output system be f(a,4, Aya? ee> Anns Sy. SQ -+e> En) The parameters E49 Sgr e+e» &p are 
the "true" values of the solution, and we wish to find an interval for each element of the solu- 
tion vector which includes the corresponding é with the desired probability. The estimates x 
are given functions of the ais, or x= &(a41, Bigs ann: The density of the 's is then 
h(x, €). If the distributions of the X's are known, confidence limits are given by 


K+ ks<é < x+ ks 


for each x, where s is the standard deviation of x and ky and ky are the critical values found 
from the distribution of x corresponding to the desired level of confidence. 

It is clear that if we desire to establish confidence intervals for the solution, the dis- 
tribution of the estimate of the solution must be known. In certain special cases which are not 
likely to be realized in practice, the form of the distribution of the solution can be known a 
priori. For example, assume that the Leontief matrix has the form 
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where the elements on the right-hand side are submatrices. The inverse of this matrix is 
given by 


— 








ie ie al eae a 
L -Ly P, Ly eee -L,, P, Ly 
-1 
0 L; 0 
0 0 i 
provided that the submatrices L,> 9 eee » Ly are square. If Lo is m x m and contains no 


errors and if errors in the remaining submatrices are independently distributed, the first m 
outputs Xj,+++,X,, are represented as sums of independently distributed random variables. 
Under certain simple conditions the distributions of the first m x's are asymptotically normal 
by the central-limit theorem [16]. The approximation will be better the larger the number of 
submatrices of type L into which the original matrix can be partitioned. 

Occasionally it may be reasonable to assume that the solution of a Leontief system is 
distributed according to the logarithmic normal distribution. This conjecture was confirmed 
by some sampling experiments with random Leontief matrices. In a case of symmetrically 
distributed errors, the normal distribution provided a good approximation to the distribution 
of the elements of the solution [17]. 


6. A NUMERICAL APPLICATION 
Assume that the true coefficient matrix is 


0.3 0.2 


0.4 0.5 


and the bill of goods is y, = yp =1. 
The Leontief matrix is 
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The determinant |I - A| = 0.27 and 


2.59259 


(I-A)? = 


1.48148 2.59259 4.07407 


1.85185 aries x, 


= 


The error matrix H is assumed to have the following properties: errors are independently 
distributed and a particular error hj; is either +0.01 or -0.01 with equal probability. The error 
distributions are symmetric with zero mean and variance 0.0001. By the use of approximation 
(4.14), 


e 1.85328 0.74224 
E[(I -A- H) ] = ? 
1.48366 2.59524 


and therefore 


E(x,) = 2.59522 , 
E(x») = 4.07890 . 
Applying (4.15) to (4.20), we obtain 
var (b, ,) = 0.002233 , 
10 var (by9) = 0.002892 , 
m 
_ var (bp) = 0.005013 , 
a. var (by) = 0.006406 , 
re) 
covar (by, > Py) = 0.002073 , 
1 is 
- covar (by » Dao) = 0.004208 . 
4 
- Applying, finally, (3.2), we obtain 


var (x,) 0.009271 , 
var (xq) = 0.019835 , 


0.09624 , 


7 (x,) 
8 (Xq) 


The expectation of the solution and the variances have also been calculated exactly by applying 
to the 16 possible images of the inverse and the solution the definitions of the mean and the 
variance. The relevant moments are 


0.12508 . 





E(x,) = 2.58581, E (x9) = 4.06413, var (x,) = 0.00794, var (X) = 0.02142. 
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One may conclude that the approximations are reasonably good. But the approximations were 
calculated on the basis of the true inverse elements. In order to avoid using the true inverse, 
fifty random matrices were obtained by shocking each element of the true matrix either posi- 
tively or negatively by 0.01 with equal probability. The solution was obtained for each of the 
50 systems of simultaneous equations, and the means and variances of the elements of the solu- 
tion were calculated. (Column 4 in Table 1.) In addition, a mean Leontief matrix was obtained 
from the sample of 50 matrices and its inverse calculated. The approximations to the expecta- 
tion and the variance also were calculated by using sample values rather than true values of the 
inverse elements (Column 3 of Table 1). The approximations are reasonable as far as the 


TABLE 1 





2 3 


Approximations 
Based on Sample | Sample Values 
Inverse Elements 


Approximations 
Based on True 
Inverse Elements 


Exact Values 








E (x,) 2.58581 2.59552 2.60008 2.59931 
E (9) 4.06413 4.07890 4.10821 4.10709 
var (x,) 0.00794 0.00927 0.00943 0.00728 
var (X») 0.02142 0.01984 0.02130 0.01560 
s (x,) 0.08911 0.09628 0.09711 0.08532 
8 (X») 0.14636 0.14085 0.14595 0.12490 























expectation of the solution is concerned. The approximations to the variances are less good 
but still reasonable in view of the fact that the percentage errors in the corresponding stand- 
ard deviations are considerably smaller. 

The difficulties of such an approach may appear to be quite formidable. In general, we 
do not possess an unbiased estimate of the true solution. The approximations to the expectation 
of the solution (which itself is not an unbiased estimate of the true solution) are based on biased 
estimates of the inverse elements. The presence of bias need not be serious if its numerical 
magnitude is small [18]. This is borne out to some extent by the hypothetical example consid- 
ered above. Further investigation into the problem may reveal more about the magnitude of 
bias. It is plausible, however, to assume that the estimates for the solution are consistent, i.e., 
as we increase the size of the sample on which the average coefficient matrix is based, X——> ¢ 
with probability one. Previously we have considered the alternative of calculating the mean 
solution rather than the solution corresponding to the mean Leontief matrix. This method can 
at least guarantee that estimates of the solution are unbiased estimates of its mathematical 
expectation. Instead of calculating a mean Leontief matrix on the basis of a sample of coef- 
ficients, we take repeated samples, obtain several sample Leontief matrices, invert each of 
them, and thus calculate the mean solution. The procedure is certain to eliminate some bias 
but is, of course, more costly in terms of computation. 

The foregoing discussion has made no use of any specific assumption concerning the 
nature of the distribution of the errors. The problem of establishing confidence intervals is 
thus quite difficult and can be solved only on the basis of empirical knowledge concerning the 
distribution of the solution. As alluded to before, the distribution of the solution might be 












































as 
tic 
the 
sti 
ar 
Cc 
er 
lis 
ap 
tic 
the 
is 


er; 










re 
3e, 
i- 


olu- 
ned 

cta- 
f the 









PROBABILISTIC ERRORS IN THE LEONTIEF SYSTEM 

















cause the scale of the diagram is too small 38 
to show the difference. B shows confidence 
intervals on the assumption of normally dis- 
tributed errors with zero mean and variance 3.6+ 
0.0001, i.e., the normally distributed errors 
are assumed tohave the same meanand var- 
iance as errors with arbitrary distributions. 
C represents confidence intervals for the Figure 1 

solution obtained by using the estimates for 

the standard deviation developed in Sec. 4 (approximation based on the true inverse and employ- 
ing the critical values of the normal distribution). D gives confidence limits by Chebyshev's 
inequality. In all cases the degree of confidence is 95 percent. It is interesting that the ellip- 
tical region B, which was obtained on the basis of the considerably restrictive assumption of 
normally distributed errors, does not give a confidence region which is uniformly narrower 
than C. As expected, Chebyshev's inequality results in wide intervals. 


approximated closely by the normal distribu- Xo 

tion if the distribution of the original coef- 

ficients is symmetric. If it is assumed that “— 

the original coefficients themselves have the 4.5+ 

normal distribution, the problem of construct- oak: 

ing a confidence region for the solution can 

be solved neatly [19]. Figure 1 shows con- 4.3+ 

fidence regions for the hypothetical system i 

under consideration, on the basis of varying 

assumptions. Point A represents the ap- 4.1r 

proximation tothe expected solution (arrived ane " . 
at by using the values of the true inverse). 

It could also represent the true solution, be- oer LZ 


3.7r 
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7. CONCLUSION 

There are some structural and sampling rationales for considering input coefficients 
as random drawings from some universe of coefficients. Under these circumstances the solu- 
tion vector of a Leontief system possesses stochastic properties. The standard deviation of 
the solution can be approximated with high degree of accuracy if (a) unbiased estimates of the 
standard deviations of the input coefficients are available and (b) the errors in the coefficients 
are sufficiently small such that the infinite series in (3.6) and (3.14) converge fairly rapidly. 
Confidence intervals for the mean solution can be established on the basis of the known prop- 
erties of the sampling distribution of the mean. A confidence region C for the solution is estab- 
lished on the basis of the conjecture that the elements of the solution vector are distributed 
approximately as the logarithmic normal distribution (which converges to the normal distribu- 
tion as the skewness approaches zero). A confidence region B is established on the assumption 
that the errors themselves are normally distributed. The intersection of the regions B and C 
is properly contained in each. 

The computational difficulties which arise in this approach promise to be quite consid- 
erable in the case of any system of reasonable dimensions. The significant methodological 
change is that the gross outputs required to satisfy a given bill of goods are not asserted as 





point estimates but are considered to be included in stated intervals with arbitrary degree of confidence. 
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ON A CERTAIN PROBLEM IN LINEAR PROGRAMMING* 


R. J. Taylorfand S. P. Thompson 
Naval Research Laboratory 


1. INTRODUCTION 

In this paper we discuss a theorem in linear programming which is fundamental to 
recent advances in design techniques for certain military vehicles. While several examples of 
such applications of the theorem appear in the classified literature, we have been unable to 
paraphrase these into a form suitable to motivate a discussion for presentation in a journal of 
this kind. At the same time we consider it important to have the purely mathematical side of 
the work published in the open literature, where it can be freely discussed and where it can 
invite unforeseen applications. 

Consider the functions 


*- 

(1.1) G= )' (i) a, (i) 
j=1 i=1 os 

and 
N 

(1.2) w=)iwed ai, 
j=1 i= 

where, for all integral i and j in the intervals 

(1.3) l<ic<a 

and 

(1.4) 1<j<N, 


x; (i) and Wj are given positive constants and where there apply the restrictions 


(1.5) 0 < aj(m) < aj(mj - 1) < ... < aj) < 1. 





*Manuscript received November 16, 1956. 
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We state and prove a theorem which determines the values of aj (i) that maximize G for 
W, subject to 


N 
(1.6) 0<W< Pa ny; Wj 
j=1 


but otherwise arbitrary. 

The linear form (1.1) measures the quantity which it is desired to maximize through 
proper choices of the variable aj (i). In the original setting of the problem, the aj (i) in turn 
measure the amount and manner of use of a particular engineering material. The x, (i) are 
constants and represent the increments in G which result from unit increments in the aj (i); 
in the original setting they are determined by unalterable features of the vehicle being designed 
and by the environment in which it is to operate. In the original setting the linear form (1.2) 
measures a penalty imposed by utilizing the aforementioned material in the manner specified 
by the a; (i), a unit increment in one of these variables imposing an increment in penalty equal 
to w.. The ordering relation (1.5) between the variables is a special feature, and insures in 
the original setting that the utilization of the engineering material is physically realizable. 

It is helpful to remark that the cumbersomeness of much of what follows stems from the 
possibility that N > 1 in (1.4). The special case N= 1 is very much simpler in all respects, 
and the reader will find that reduction to this case will frequently clarify many of the com- 
plexities which must be considered. 


2. DEFINITIONS 

The following definitions divide the interval (1.6) for W into nonoverlapping subintervals 
and allow the theorem to be stated in the form most easily recognized by those having knowledge 
of the applications. Other readers may prefer to read first the statement of the problem and the 
theorem at the beginning of Section 4, which result from a formal change of variables. These 
are more compact than the corresponding forms which we are now presenting but, from the 
point of view of the existing applications, also somewhat artificial. 

Define for each integer j in the interval (1.4) and for the integers a ands which are 
subject to 


(2.1) 
but which are otherwise arbitrary, the quantities 


a 
(2.2) Xj(ap)= D> x (i), 
i=B+1 


X; (a, 8) 
C;(a,6) pe havens ° 


bed (a-B) 


It may be seen that (2.2) represents the increment in G, which results from increasing 
aj (i) from zero to unity for all i in the interval 8 + 1< i < a, and that (2.3) is the ratio of 
this increment in G to the corresponding increment in W. In the original setting of the 
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problem, Cj (a,8) may be interpreted as a measure of the efficiency which accompanies use in 
the foregoing manner of the particular engineering material with which that setting deals. But, 
regardless of the setting of the problem, these quantities play a central role in the theorem, 
and the purpose of the remainder of the Section is merely to order them in a way which is con- 
venient for the statement and proof. 

Associate with each value of j the sequence of successive positive integers: 


(2.4) 8; = a, By o0es 


terminating with an integer Si >. 1, whose definition will follow. With each value of 8; associate 
the integer 


(2.5) k,(s;) > k,(0) = 0, 


which is defined by requiring for arbitrary integral i in the interval 


(2.6) k; (8; - 1) <i< nj 
that 
(2.7) Cc [k; (8), k; (8; -1)] > Ci [i, k; (8; -1)}. 


Finally, define Sj as the integer such that 


(2.8) k; (8;) = nj. 

It will be seen that for each sj 

2. .(s;- .(s. .(S:) = n,. 

(2.9) kj (sj-1) + 1 < (8) < k(S)) = 0; 
Irrespective of the values of j, order the totality of quantities 

2. .[k. (s; .(s;- . ‘ j 

(2.10) Ci[kj(s), k(sj-1)], 1 <8; <8, 1<i<N 


into a sequence of nonincreasing magnitudes and place the quantities ordered in this sequence 
into correspondence, respectively, with the integers 


Mz 


(2.11) 340% 


S;. 
j J 


1 


By means of this association there corresponds to each integer q in (2.11) a value of j, to be 
denoted by j(q), and a value of s,, to be denoted by 


For each j in (1.4) and for each q in (2.11), define v5 (q) by 
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85 (q) (®, if j S j(q), 
9 ;(a) = 
0, if j 4 ia). 


Then for each j in (1.4) and each q in (2.11), define by P; (q) the largest member of the 
set 


(2.14) 0;(a) , 


where 


(2.15) l<ac<q, 
i.e., the member such that for every @ in this interval 
(2.16) p;(a) > oj (a). 


From Appendix I and the ordering specified in connection with Eq. (2.10); it will be 
seen that 


(2.17) 85 (q+ 1) (4 +1) = Pi (q+ 1) @ +1 


and therefore that 
(2.18) Pj (q+1) (4+ 1) = Pj (q+1)(@ + 1. 


For each q in (2.11), we finally define W(q) by 


N 
(2.19) wea) = 2, wk a), 
j= 


which values divide the interval (1.6) of W into nonoverlapping subintervals. 
The following theorem applies to the subinterval 


(2.20) W(a9) < W < W(ap+ 1), 
wherein it is to be understood that Io is initially selected arbitrarily from (2.11) but is there- 


after held fixed. This particularization allows a considerable simplification of notation. For 
the remainder of the paper we set 


(2.21) res (dp + 1) 
and for all j in (1.4) shall understand the argument of Pj to be Qo - In addition, we shall hence- 


forth omit the subscript j on k; whenever it is evident from the argument. 
With these definitions and recalling the ordering specified in (2.10), we can say 
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(2.22) C[k@), kp; - 1] > C,[k@,+ 1), kP,)] > Cli, kD], 
for arbitrary j in (1.3) and for arbitrary i in the interval 

k(p;) <i< ~ ’ 
an inequality we will find necessary in the proof of the theorem. 


3. THEOREM 
For W arbitrarily selected from the interval 


W(d9) < W < W(qa9+ 1), 


the maximum value of G is given by 





Gt) max G90] = ep et) Rp] tReet) Be + X x; [k(p), 0], 


J 
and the maximizing aj (i) are given by: 


cr 








W - 2 wy k(p)) 
a_(i) = =e ae * k(p,) < i < k(p,+1), 
(3.2) a,. (i) = 0, for k(p,.+ 1)<i< n,; 
aj (i) = 1, for 1<i<k(p)), and all j, 
| yM=0, for k(p;)<i<nj;, andalljyr. 


Remark: Regarding the asserted maximizing aj (i)'s as functions of W, the reader may see 
from (3.2) and (2.14)-(2.19) that the theorem implies for 


N 
(3.3 n. w; > Wy > W, > 0 
) Ps ie 2 1 
and that 
(3.4) (a; Ow, = (a; Ow, 


for all i in (1.3) and all j in (1.4). 
The foregoing statement of the theorem insures that the assertions contained within it 
are made precisely. Unfortunately, the same form also conceals the major ideas in these 


assertions, which are quite simple. Therefore, at the risk of being imprecise, we describe in 


words what the optimum assignment of aj (i) is asserted to be. 
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Compute all of the quantities defined in (2.3). Consider the largest of these and ask 
whether or not the associated a;(i) can all be set equal to unity, all other a(i) being zero, 
without violating condition (1.5). If the answer is no, consider similarly the next largest value 
of C.(a,8); if yes, assign the associated a,(i) the value unity and maintain these assignments 
for the rest of the process. From the list of C; (a, 8) strike the one whose a.(i) have just 
been assigned. Repeat the foregoing operations, beginning with the largest C;(a,8) which is 
still listed, remembering that certain of the a, (i) have now the value unity. As one continues 
in this way, many of the a. (i) will be assigned the value unity; these appear above as the third 
member of (3.2). A stage in the process will finally be reached where unit values for the 
relevant variables are permitted by (1.5) but where (1.2) would cause the preassigned value of 
W to be exceeded. When this occurs, set these variables equal to one another and to a value 
sufficient to satisfy (1.2); these appear as the first member of (3.2). The process of assign- 
ment is now complete, and the remaining, zero valued, variables in (3.2) correspond to C; (a, 8), 
which are sufficiently small that the process is completed before they are reached. 


4. PROOF 


In order to simplify the proof of this theorem, we introduce the following transformation 
of variables. Define for all j in (1.4): 


b; (0) = 1 - aj (1) ? 





| By mp = aj (n;) : 


Inequalities (1.5) are satisfied for arbitrary aj (i) if, and only if, 
n 


j 
dy bj (i) = 1, for arbitrary j, 


and 
(4.3) b; (i) >.0, for arbitrary j andi. 


From Eqs. (4.1) through (4.3) it may be seen that the foregoing transformation possesses 
a unique inverse, namely, 


aj 


(4.4) a; (i) = 2d, by(y), l<cicnm, 1<icN. 


It may be seen that under this transformation G, Eq. (1.1), and W, Eq. (1.2), become, 
respectively, 


e : 

(4.5) G(b) = b, (i) X,[i, 0 
foi fea * _ 
and 
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. 
(4.6) W = W; > ibj (i). 
j=1 ° i=1 


From the use of Eqs. (3.2) and (4.1), it will be seen that proving the theorem is equivalent to 


proving that G(b) is maximized by the following values for b (i): 


W - Dw; k(p.) 
j J J 


w,.[k(p,.+1) - k(p,)] ’ 





(4.7) b..[k(p,.+1)] = 


Pa Wj k(p;) . w,k(p,.+ 1)-W 
(4.8) bt. [k(p,)] = 





w,[k(p,+1) - k(p,)] 

(4.9) b).(i)= 0 for arbitrary i 4 k(p,+1), k(p,), 

(4.10) b; [(p;)] = 1 for arbitrary j¥r, 

(4.11) bj (i) = 0 for arbitrary i # k; (P;) and arbitrary j # r. 


We can now prove G(b') is the maximum G(b). In order to accomplish this, we shall 
prove, for arbitrary b, that: 


(4.12) G(b') - G(b) > 0. 


The proof from here on is entirely formal. While lengthy, it has the advantages of being 
rigorous, entirely self-contained, and of requiring only elementary algebra as a prerequisite. 
The inequalities contained in the definitions are used to derive additional inequalities, which are 
then arranged to show that the left-hand side of (4.12) can never be negative. The necessary 
expressions are cumbersome, but will be found to undergo a significant reduction in complexity 
for the special case N=1. The reader who wishes to proceed further is urged to consider this 
case by first suitably reducing the following general treatment. 


Substituting the b' values, (4.7) to (4.11), into (4.5), we obtain: 


* Wj k(p;) + w,K(P,+ 1)-W 
G(b') = 





X_|k ,0 
w,[k(p,*1) - k(@,)] ital 


W - Dw; k(p,) 
ie A e.. 





X, [k(p,+1), 0] + Pa X; [k(p,), 0]. 


+ 
w,.[k(p,+ 1) - k(p,)] 
By noting the definition (2.2), we can write 


(4.14) X,.[k(P,+1),0] = X,[k(p,+1), k(p,)] + X,[k(P,), 0]. 
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Substituting this in the second term of (4.13) and collecting terms, we then obtain 


w-2 w; k(p)) 


j 
(4.15) G(b') = v,[k@,* 1) - k,] X, [k(p,+ 1), k(p,)] + . X; [k(p;), 0}. 





If we consider G(b) for general bj (i) as presented in Appendix II, we can write 


aj 


N 
coo eam ® 7) 
j=1 


i-k,(1}+1 


k;(1) 
-W; a ib;(i) C,[i,0]+ C;[k,(2), k;(1)] w; a -k,(1) 


7 
- D* faye) - 0) no 
i=k;( }1 


(2) 

-w; [i - K(1)] b\(i) Ci, Ky(A)]+ C,L,(3) - 2], ee - k,(2) 
i=k,(1}+1 

n 


a * 


j 
[,(3) - k,(2)] 0) 
i=k,(3) +1 


k,(3) 
” bad je [i-k,(2)] b,(i) C;[i,k\(2)]+---+ C;[K(p;), k(p;-1)] a op ws k(p;-1) 


i=k,(2}+1 


. 


- D*  [e(p,) - k(p;- 0] “0 
i=k(p;}+1 
kip;) » | 
-wi »: [1-k(p,-1)]b,@)C,[4,(9,-1)] -w, Zz [i-K(p,)]b,(i)C; [i, tp) 
isk(p;-1}+1 isk(p;}+1 
+ C,[k(p,.+1), k(p,)] [w - . wk() | . 


This equation can be written as follows: 
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N . 


Pj J 
G(b')- G(b) = >" m Is [k(s,),k(s,-1)] w, a K(s;-1)- | mapae-IH(0 


isk(s;}+1 


k(s;) 
“7 2: [i-k(s;-1)]b,(i) cf 
isk(sj-1}+1 


[i-k(p;)]b;(i) C,[i, K(p;)] > + C,[k(p,+1),k(p,)] lw-2 L wik(p; |: 


From (I.6) and (2.7) we can write 


k(s;) k(s;) 
(4.18) - ba Zz. [i-k(s ;-1) ]b;(i) C;li,k(s;-1)] = bal ¥. [i-k(s;-1)] b (i) C;[k(s;),k(s;-1)], 


isk(sj-1}+1 isk(s;-1}+1 
and from (1.6) and (2.22), 


n,. Nn, 


J 
(4.19) - wy [i-k(p,)] b,(i) CyLi,k(p,)] > -w) )* — [i-k(p,)] bj(4) C; [k(p,+1), k(P,)] - 


i= =k(p;}+1 


Nn: 


j J 
G(b') - G(b) 2. >: C;[k(s;),k(s;-1)] a {Mapa >: [k(s;) - k(s;-1)] bi(i) 
i=l s.=l i=k(S;}+1 


txt} + C,[k(p,+1), k(p,)] 


isk(Sj-1)+1 


7 
\w Yayo) Dy pm tatoo | 


j=1 isk(p;}+1 


§ The coefficient of C. jlk(s, ), k(s; -1)] is positive (see Appendix III), hence we can apply (1.6) to 
these coefficients. ieusine in nous the inequality (1.10), we can write 
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N » j 

G(b")- G(b) > ) | Clk), tp,-1)] wD 4 R(s;)- k(s,-1) - [e(s})-k(s;-1)] 0,0) 
j=1 s.=1 i : 

J 


k(s;) 


N 
[i-k(s,-1)]b,(i) + C_[k(p ot e9{ “fe w; k(p,) 


i=k(s;-1}+1 


N > 
Dr DD. f-k@p)1 bw 
i-k{p,)+1 


Note now that the coefficient of C,[k(p.), k(p,-1)] in (4.21) is positive for each j, being 
the sum of coefficients of C;[k(s)), k(s;-1)] in (4.20), each proved to be positive in Appendix if]. 
This coefficient can be written as 
k(p;) - 

(4.22) kp) - D) ibf- D> km) di), 
i=1 i=k(p;}+1 

from Appendix IV (IV.3). Again applying (1.6) and (2.22), we can write 


k(p;) 
G(b') - G(b) > C,[k(p,+1), k@,)] y k(P,) w; - ye > ib,(i) 


(4.24) 
=0, by (4.6). 


Hence, G(b') - G(b) >.0, so the b''s defined by (4.7) to (4.11) result ina maximum G(b). By 
use of the transformations given by (4.1), the maximizing a's can be found from the b' 's and are 
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N 
w* Z wjk(p;) 
j=1 

a(i) = for k(p,)+1< i< k(p,+1), 


w,[k(p,+1) - k(p,)] 





a_(i) = 0 for k(p.+ 1) <i <n,, 


a; (i) = 1 for 1<i<k(p)), and all j, 


a; (i) = 0 for k(p;) <i<nj, andall jfr. 


The theorem is then proved, and we note that the form of the maximum G, G(b'), given 
by (4.15), is the one given by the statement of the theorem (3.1). 
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APPENDIX I 


From (2.7) and (2.3) we can write 


X; [k(s;), k(s;- 1)] é X; [i - k(s; -1)] 
k(s;) - k(s; - ” ~ ts k(s; -1) 





(1.1) 
for arbitrary i in the interval k(s;-1) <i< nj. Hence, 


X; [k(s;), k(s; -1)] Xj [k(s;+ 1), k(s; -1)] 
k(s;) - k(s;-1) > k(s,+1) - k(s, -1) 








(1.2) 


X; [k(s;+1), k(s;-1)] [ k(s;+1) - k(s)) | 





k(s;+ 1) - k(s;) | k(s;+ 1) - k(s;~- 1) 








x; [k(s;+ 1), k(s;-1)] [ k(s;+ 1) - k(s;-1) + k(s;-1) - k(s;) 
k(s;+ 1) - k(s;) q k(s;+ 1) - k(s; -1) 


(1.5) 








Xj [k(sj+1), k(s;-1)]  [ k(s;) - k(s;-1)] X;[k(s;+1), k(s;-1)] 
~  k(8;+1)- (8) | (#1) - K(s)) | k(s;+1) - k(8;-1) 


Note that if 
23 
and K is an arbitrary positive constant, then 
KZ, < KZ, 
and 
(1.6) - KZ, >-KZ,. 


From this and (1.2), 








- et - | X; [k(s;+1), k(s; -1)] — - =tes| X;[k(s;), k(s; - 1)] 
I.7) - > a 


k(s;+ 1) - k(s;) k(s;+ 1) - k(s; - 1) k(s;+ 1) - k(s;) k(s;) - k(s;- 1) 


Then, if we recall (2.2), (1.5) can be stated as follows: 
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(1.8) 





X;[M(s;), k(s; -1)] , X;[k(s 5+ 1), k(s;)] X;[k(s;), k(s; -1)] 7 k(s;) - k(s; -1 X;[k(s;),k(s;-1)] 
K{(s;)-k(s;-1) ~~ k(s;+1)-K(s,) r K{s,+1)-K(s,) ) 


k(s 5+ 1)- K(s; k(s;) - k(s; -1) 
7 X; [k(s;+ 1), k(s;)] 


ated k(s;+ 1) - k(s;) 





Hence, given (I.1): 


(1.10) C; [k(s;), k(s; -1)] >. C; [k(s;+ 1), k(s;)] . 
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APPENDIX II 


Consider the general term in G(b), Eq. (4.5), for arbitrary b and j. It can be expanded 


as follows: (11.3) 
n, K,(1) (2) - 
b oe b; (i) X[i,0] = >: b; (i) X; [1,0] + b; (i) X; [1,0] +. . . 
i=1 i=1 i=kj(1}+1 
k(p;) * — 
+ >: b; (i) X; [1,0] + 2. b,(i) X; [3,0] - 
i=k(p;-1}+1 i=k(pj}+1 r 
From the definition (2.2), we can write ) 
l- 
(1.2) X, [1,0] = X,{k,(1),0] + X; [kj (2), k, (1)] +... + X;[k(s,), k(s;-1)] + Xfi, k(s)], 
for i> k(s;), but otherwise arbitrary. 
Using (11.2), we can further expand (II.1) as follows: 
a 
b.(i) X.[i,0 
X Pe 
(1.4) 


— sma 3 x,[k(1),0) 5) %,[i,«(0)] 


= ) Ibi) + >: KD) ay >: [i-K,(1)]b,(i) iq) 
i=i i-ky(1}+1 isk (T+ 


k;(3) X;[k,(1),0] k;(3) X;[k,(2),k,(1)) 
Pw 2, Uig(2)-Ky(a)]by(0 K@)-j) 


i-k,(2)+1 i=kj(@)}+1 


k,(3) K,[isk(2)] Mp) X,[k,(1),0] 


+ D) [K@b) tet DCO 


i=k,(@)+1 i-k,(2) i=k(p;-1)+1 K(1) 


xt xen) “PP — [K(3),42)] 
+ [k,(2)-k,(1) (i) Kar KI) * alleles [,(3)-K,(2)]b,(i) “Ke 
=K(p 





isk(p;-1}+1 


k(p.) n. 
j X.[i,k(p,-1)] ) X,[k.(1),0] 
ot bs [i-k(p,-1)],() | 7 is Hono 
i=k(p;-1}+1 J i=k(p;}+1 
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Nn: Nn. 
X;[k,(2),k;(1)] . X;[k,(3),k;(2)] 
+ a [k,(2)-k.(1)] b.(i) : + [k.(3)-k.(2)] b.(i) 
i=k(p;}+1 wr ' k;(2)-k;(1) isk(p;)}+1 _ ' k;(3)-k;(2) 








Nn: 


X,[k(p,),k(p,-1)] X [Mp] 
i 


" 
tenet [k(p.)-k(p.-1)] b.(i) + [i-k(p.)]b.(i) ——_—_.. 





By collecting terms, we then obtain 


ny k;(1) 


n. 
X [1,0] X,[k,(1),0] 
Z b,(i) X,[i,0] = } i by(i) * + y Kono a 
i=l i=1 i=k,(1}+1 
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The coefficient of C;[K(s,), k(s;-1)] in Eq. (4.20) is positive. This coefficient can be 
written by using (4.2) as: 


ni K(s)) 
[k(s;)-k(s;-1)] yD b;(i)-[k(s;)-k(s;-1)] bj(i) = 2 [i-k(s;-1)] bi(i) 
i=0 i= s;}+1 i=k(s;-1}+1 


k(s;) k(s;) 


= [k(s,)-m(s;-1)] )) BYi)- — D?—_[i-Wfj-1)] b(t) 
i=0 i=k(S;-1)}+1 


= k(s;) Ly b,(i)- k(s,-1) i. b,(i)+ k(s,-1) be b,(i) - ta ib,(i) 


i=k(s;~1)+1 i=k(s;-1)}+1 


5 k(s,-1) k(s.) 
= (s,) bin b,(i)+(s,) b(i)-k(s-1) D' vfi)- Dib, 


ick; bias i=0 i=k(s;-1}+1 ; 


k(s 5-1) k(s 5) 
= [k(s,) - i(s,-1)] he bi) + >: Ds) - 4] B,C) 
i= i=k(s;-1)+1 


>.0 from (2.9) and (4.3). 





ON A CERTAIN PROBLEM IN LINEAR PROGRAMMING 


APPENDIX IV 
Coefficient of C;[k(p)), k(p; -1)] in Eq. (4.21) 


From (4.21), this coefficient can be written 
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NOTES 
PROBLEMS 
It has been suggested that this journal might serve as a medium for publishing problems in 
the area of Logistics — with the idea that other persons might be interested in those problems 
and might submit comments. Readers are invited to submit brief statements on applied and 
theoretical problems in Logistics. Address letters to Managing Editor, Naval Research Logistics 
Quarterly, Office of Naval Research, Washington 25, D. C. 


The editorial note for the paper “Dynamic Problems in the Theory of The Firm,” by 
H. M. Wagner and T. M. Whitin, which appeared in Volume 5, No. 1, (March 1958) should read 
as follows: 

This article was published as Appendix 6 in the second edition (1957) of the book entitled 
“The Theory of Inventory Management,” by Thomson M. Whitin. It is being published in this 
journal in order to reach an audience which might otherwise not be reached in the field of 
logistics. The Princeton University Press, who publish this book and who are the copyright 
owners, have granted permission to reprint the above material in the Naval Research Logistics 
Quarterly. 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Editor items of general interest in the 
field of logistics 


DR. I. STAKGOLD has relieved Dr. F. D. Rigby as Head, Logistics Branch, Office of Naval 
Research. Dr. Rigby is Director, Mathematical Sciences Division, Office of Naval Research. 
Dr. I. Stakgold has joined the board of editors of the Naval Research Logistics Quarterly. 
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RECENT PUBLICATIONS 


EDITOR’S NOTE: Listing of a publication in this section is for record and reference only and 
does not constitute an endorsement of point of view or advocacy of use. 


U.S. ARMY IN WWII. SPECIAL STUDIES - REARMING THE FRENCH. By Marcel 
Vigneras. Published by the Office of Chief of Military History, Department of the Army, 
March 1958. On sale by Superintendent of Documents, Washington 25, D. C. $4.25. 


When French national forces went back into the war against the Axis powers in mid- 
November 1942, their re-emergence on the battle fields of Europe after a temporary eclipse 
heralded the start of an early large-scale American experiment in mutual aid as well as the 
resurgence of a gallant ally. Since 1940 French participation in the common struggle had been 
limited to General de Gaulle's small band of Free French, armed and maintained by the British 
and stubbornly fighting on their side. Marcel Vigneras' Rearming the French tells the story of 
how a much larger French force, under the leadership, sometimes challenged, of the redoubt- 
able General Henri Giraud, swelled the ranks of Allied military manpower, helped tip the bal- 
ance of power in the Mediterranean, and ultimately assisted in the liberation of Europe. 





A French combat force that eventually numbered nearly 400,000 men was equipped, 
trained, and put into action on land, in the air, and on the sea—all in record time. Obviously, 


an undertaking of such magnitude was complicated by innumerable obstacles, and Dr. Vigneras 
does not hesitate to point them out in the course of his analysis of how the venture fared. Bas- 
ing his account on a wealth of American and French documentary sources bolstered by many 
personal interviews with and letters from participants, the author ranges over a wide series of 
interacting influences on the course of French rearmament. Among the significant factors dis- 
cussed are the Giraud-de Gaulle tug-of-war and other French political developments; the posi- 
tions of the Americans, the British, and the French themselves regarding the undertaking; the 


competing global demands for shipping and armament; and the limitations of French technical 
personnel. 


Because mutual aid has continued to the present day as an integral part of United States 
policy, Rearming the French will prove a boon to those who may have to carry out similar 
wartime enterprises in the future. For it is in effect a case history of problems that are likely 
torecur. Among them Dr. Vigneras highlights the language barrier; differences in national 
customs, food habits, and clothing sizes; requirements in liaison, armament, and training per- 
sonnel; stock shortages and surpluses; control of deliveries and distribution of equipment; and, 
to complete a very slim sampling, requirements in support and service troops. 





Opening with a brief review of the assistance similarly extended by France to an un- 
prepared America in 1917 and 1918, the narrative covers the beginnings of the rearmament 
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program in North Africa and proceeds through the eventual entry of the North African forces 
into combat in the Mediterranean, Italy, and France. The latter portion of the book takes up 
the thorny question of liberated manpower in France and describes the abortive French attempt 
to set up a Far East Expeditionary Corps. Dr. Vigneras closes his work with a concise totting 
up of the debits and credits stemming from the venture. There is no question that it paid off- 
in dollars and cents, in the restoration of prestige and pride to an old and valued friend and 
ally, and in thousands of American lives. 

(Reviewed by Ruth Stout) 


THEORETICAL WELFARE ECONOMICS. By J. de V. Graaff. Cambridge University 
Press, New York, 1957. $4.00. 


This is a frustrating book to review. On the one hand, it is written by a very able 
economist, and it is undoubtedly the best available summary of welfare economics in print. 
Most of the strands in current thought are presented and the underlying assumptions brought 
out with some care. Yet there is a curious lack of real life in the book, perhaps to be explained 
by the author's low opinion of the whole subject, as brought out in the last two chapters. There 
is a strong tendency to emphasize the limitations of possible conclusions and the strength of the 
assumptions needed to arrive at them. 


Mr. Graaff accepts pretty completely the Samuelson-Bergson version of welfare 
economics, which has two parts: the Paretian marginal conditions for an optimum in the weak 
sense ofa vector maximum in the space of individual's utilities, and the social welfare function 
The logical difficulties in the latter are not faced up to. In the former, the reader is led throug 
several stages relating to technology, tastes, and potential welfare (a combination of the two 
preceding). The chief novelty, apart from the care and thoroughness of the exposition, is the 
introduction from the beginning of external effects. But in accordance with the author's general 
attitude, they generally are studied only to show that previously accepted results do not hold. 
They are introduced only in a very general way, the output of one firm being made to depend 
upon the output of another, without any more specific analysis which might supply the basis for 
a theory of achieving optimal states in the presence of external economies. 


The author makes no reference to and no use of the tools of activity analysis, linear or 
nonlinear, in presentation of technological efficiency, so the whole presentation appears some- 
what dated. (As a matter of fact, I would guess from internal evidence that the book was com- 
pleted not later than 1950.) 


The author goes on to discuss the problems raised for welfare economics by investment, 
indivisibility, and uncertainty. The issues with regard to the first are treated fairly well, but 
his insistence on a nonmarket solution depends to some extent on the weakness of his assump- 
tions. (For example, in the von Neumann model of the expanding economy, the optimal rate of 
growth can be achieved through the market, and one can easily find weaker sets of assumptions 
for which it is true.) The discussion of indivisibility rehearses the well-known considerations 
competently; that of uncertainty is useless. Since the author has rejected the von-Neumann- 
Morgenstern utility construction (without any really serious argument), he has no tools to 
analyze the problem and projects this lack of tools to welfare economics as a whole. 
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The chapter on foreign trade, which originally appeared as an article, is, on the con- 
trary, a highly original clarification of the problem of the optimum tariff and is clearly the 
high spot of the book. 


The final chapters deal with marginal cost pricing, the evaluation of social income, and 
a set of conclusions; they serve chiefly to bring out the author's distrust of the whole subject 
and emphasize a rather naive faith that positive economics untroubled by welfare considerations 
will serve to determine policy. 
(Reviewed by Kenneth J. Arrow) 


GAMES AND DECISIONS, INTRODUCTION AND CRITICAL SURVEY. By R. Duncan Luce 
and Howard Raiffa. Wiley, New York, 1957. 509 and xix pp. $8.75. 


The expressed intent of this book is "to communicate the central ideas and results of 
game theory and related decision making models unencumbered by their technical mathemati- 
cal details." In this the authors, both of whom are mathematicians, have succeeded remarkably 
well; for despite the book's sweeping scope, there is nothing superficial about its presentation. 
It is indeed a critical introduction which clearly denotes both the underlying assumptions and 
some conceptual weaknesses and unsolved problems, as well as the real contributions of game 
theory. Luce and Raiffa's explicit and discerning discussion of these objections, most of which 
are not strictly mathematical, is a very noteworthy effort in helping the social scientist to 
understand some of the problems which confront him daily in his own work. 


Besides a superb exposition of the fundamentals of game theory, as first expounded in 
John von Neumann and Oskar Morgenstern's Theory of Games and Economic Behavior, the 
authors have also incorporated into their text the essence of the more recent contributions in 
the abstract theory. They have also noted some of the contributions and applications of game 
theory to statistics, economics, operations research, and military strategy. Both these tasks 
have been accomplished not only with clarity but with a minimum loss of understanding which 
normally results from any condensation. 





This entire treatment has been skillfully reduced to 370 pages of text and a somewhat 
more mathematical appendix of 114 pages. The former section, which can be easily com- 
prehended by any interested scientist, is divided into 14 chapters. The first three chapters are 
devoted to an introduction to the subject, including both the concept and examples of conflicts 
of interests, an informal characterization of a game, the difference between a game's extensive 
and normal forms, and the idea of a pure strategy. This section also includes a slight digres- 
sion which outlines linear programming as an example of decision-making under certainty and 
notes the real assumptions and commonly held fallacies of modern utility theory. 


Chapters 4 to 6 concern the two-person game. The first of these chapters contains a 
lucid exposition of the basic principles of game theory and covers such considerations as 
strictly competitive games, equilibrium pairs, the minimax theory, and mixed strategies for 
zero-sum games. The next chapter adapts these same concepts to the nonzero-sum game and 
illustrates these ideas with several interesting examples. It also introduces the idea of a 
solution for a noncooperative game. In the following chapter, which is devoted to the two-person 
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cooperative game, the solution takes on added significance, since varied restrictions on pre- 
play communication, etc., destroy the uniqueness of the solution. Consequently three versions 
of a "solution" are discussed at some length. These are those proposed by von Neumann and 
Morgenstern, by Nash, and by Shapley. 


Chapters 7 to 12 are devoted to n-person game theory. Because of its inherent com- 
plexities, this area has not as yet been as rigorously defined as that of two-person theory, but 
Luce and Raiffa present with remarkable clarity the progress in this sector. In Chapter 7 they 
discuss the intriguing ideas of behavioral and composite strategies, as well as the more formal 
concepts of von Neumann-Morgenstern. Chapter 8 introduces the idea of side payments and is 
largely concerned with strategic equivalence, the characteristic function, imputations, and the 
core. The complexity of the n-person theory is demonstrated by Luce and Raiffa, who have 
devoted three chapters (9-11) to "solutions" and related ideas, and again in chapter 12 they 
note the obstacles encountered in trying to apply abstract n-person theory to real game 
situations. 


Chapters 13 and 14 survey two important decision-making problems. The first of these, 
individual decision-making under uncertainty, is very closely related to game theory. As 
pointed out by Luce and Raiffa, modern statistical decision theory is essentially a simple two- 
person game played against nature. Chapter 14 deals with a totally different problem—group 
decision-making or the Arrow problem. 


The appendix is considerably more rigorous, but even so extensive formal mathematical 
training is not a prerequisite. It covers much of the same material as the first six chapters an, 
in addition, some of the more recent advances in two-person theory. The most interesting of 
these concern games with infinite pure strategies and also an attempt toward a more dynamic 
theory by sequentially compounding two-person games. 


This book is not only an excellent introduction for the layman who may not wish to read 
the appendix but is also a valuable tool for the expert who surely will desire to do so. For 
either reader it should prove to be both informative and enjoyable. 

(Reviewed by F. R. Shupp) 


QUEUES, INVENTORIES, AND MAINTENANCE-—THE ANALYSIS OF OPERATIONAL 
SYSTEMS WITH VARIABLE DEMAND AND SUPPLY (Publication No. 1, Operations Research 
Society of America). By Philip M. Morse. Wiley, New York, 1958. 202 and iv pp. $6.50. 


Queues, Inventories, and Maintenance is an elementary exposition of queueing theory 





and its applications to waiting-line problems in operational systems. By means of a series of 
concrete problems, it introduces the main concepts and definitions of the field. These prob- 
lems should also supply the reader with a real "feel" for various phenomena and general char- 
acteristics of queues. The book is particularly timely in bringing an important aspect of re- 
search in operational systems well within the reach of an undergraduate mathematical back- 
ground. Its subject matter should set a minimum, but effective, standard for the growing 
discipline of operations research in this area. 
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It seems very appropriate that Professor Morse author this first book of the ORSA 
series. Aside from his personal influence and contributions to operations research, his book 
exemplifies much of what seems worthwhile in this field. First, the subject matter itself has 
emerged as a primary tool and promises to be an opening wedge into more sophisticated treat- 
ments of operational systems as stochastic processes. Second, Professor Morse furnishes a 
striking illustration of an applied scientist''at work. The bookis a worthwhile study in getting a lot 
of mileage out of modest equipment, through fluency in model building which exploits readily avail- 
able tools of analysis, and imaginative use of such models in economic and systems analyses. 


In spite of the excellence of the exposition and Professor Morse's delightful style, there 
are two quite serious shortcomings to the book. While there is no pretense of developing 
theoretical results, the reader is required to use such results simply in order to work prob- 
lems. Many key theoretical results which Professor Morse uses are not stated explicitly; and 
an uncritical reader may not even realize that the solution to a problem hinges on them. For 
example, given P = (P;;), a (stationary) Markoff process with transition probabilities, Pi; from 
state i to state j (of a denumerable set of states), one can see clearly that a stationary prob- 
ability (frequency) distribution of states, given by x = (x;), will satisfy the relations 

Xj = 2 x Piy 2X =1, x, > 0. 
But with an arbitrary process P, it is not obvious that: 
(a) The relations of (1) have a solution at all 
(b) If the relations of (1) have a solution, that this has any meaning for the 
process—e.g., that the process converges, in some sense, from an initial 
state to such a distribution of states. 
Theoretical results and conditions which allow one to solve a set of relations such as (1) and 
interpret their solutions as stationary distributions of states are central to many of the prob- 
lems dealt with, yet they are not referred to in the book. 


The other major shortcoming is that little or no attempt is made to integrate and ref- 
erence the book and its contents into the substantial literature of mathematics and statistics 
that surrounds it. Feller's "Probability Theory and its Applications" is an admirable refer- 
ence for the theoretical bases of the book (it is listed in the bibliography). But the reader often 
has to recognize an implicitly stated theorem (if he can) and then find it in the literature with 
n0 help from the author. A minor point is in the cavalier treatment of certain mathematical 
quantities—such as omitting a remainder in a Taylor expansion. Professor Morse does not 
get into trouble thereby, but less sophisticated readers may, indeed. 


These objections seem worth the noting because Queues, Inventories, and Maintenance 
is, in its timeliness and exposition, such a good book. The reader will find the potentials and 
general approach of ‘queueing theory expounded very nicely, but, unfortunately, few of its pit- 
alls and few clues for extensions in its power are given. 

(Reviewed by Harlan Mills) 
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